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We investigate theoretically a diffusion- limited reaction between a reactant attached to a Rouse 
polymer and an external fixed reactive site in confinement. The present work completes and goes 
beyond a previous study [T. Guerin, O. Benichou and R. Voituriez, Nat. Chem., 4, 268 (2012)] that 
showed that the distribution of the polymer conformations at the very instant of reaction plays a 
key role in the reaction kinetics, and that its determination enables the inclusion of non-Markovian 
effects in the theory. Here, we describe in detail this non-Markovian theory and we compare it 
with numerical stochastic simulations and with a Markovian approach, in which the reactive con- 
formations are approximated by equilibrium ones. We establish the following new results. Our 
analysis reveals a strongly non-Markovian regime in ID, where the Markovian and non-Markovian 
dependance of the relation time on the initial distance are different. In this regime, the reactive 
conformations are so different from equilibrium conformations that the Markovian expressions of 
the reaction time can be overestimated by several orders of magnitudes for long chains. We also 
show how to derive qualitative scaling laws for the reaction time in a systematic way that takes into 
account the different behaviors of monomer motion at all time and length scales. Finally, we also 
give an analytical description of the average elongated shape of the polymer at the instant of the 
reaction and we show that its spectrum behaves a a slow power-law for large wave numbers. 

PACS numbers: 02.50.Ey,82.20.Uv,82.35.Lr 



I. INTRODUCTION 

Among transport-limited reactions, reactions involving 
polymers play an important role and have been widely 
studied, both experimentally [1-7 and theoretically [8]- 
[16]. When a reactant molecule is attached to a polymer, 
its interaction with the whole polymer chain results in 
a complex motion that can be subdiffusive [TTJ [18] and 
leads to non-trivial reaction kinetics [13] [19] [20] . Under- 
standing polymer reactions is useful for biologically rele- 
vant problems such as the kinetics of hairpin or loop for- 
mation in nucleic acids [TH3] or the folding of polypeptide 
chains [5ti7][2T|. In these examples, the monomers belong 
to the same chain. In this paper however, we focus on 
intermolecular reactions that occur between monomers 
of different chains or between a single monomer and an 
external reactive site fixed in a confining volume (Fig. [T]), 
as in the case of the search of a pore or a catalytic site in 
a confining cavity during gene delivery or viral infection 
[22] [24]. 

The theoretical description of polymer reaction kinet- 
ics in the diffusion controlled regime is complicated by 
the structural dynamics of the chain, which implies that 
the motion of a single monomer cannot be described as a 
Markov process, as the other monomers of the chain play 
the role of "hidden degrees of freedom" . For this reason, 
the determination of the reaction kinetics is a difficult 
task, even in the simplest case of a Rouse chain model 
for the polymer that is considered in this paper, where 
hydrodynamic and excluded volume interactions are ne- 
glected [25]. The first theoretical approaches of poly- 
mer reaction kinetics overpassed this difficulty by doing 
Markovian approximations, either by replacing the whole 




FIG. 1. (color online) Sketch of the problem investigated in 
this paper. A reactive monomer is attached to a polymer in 
position p in the chain of N monomers. We calculate the 
mean time T for the reactive monomer to reach a reactive 
region of radius a in a confining volume V for the first time, 
given that the initial distance between the reactants is Ro, 
and that the internal degrees of freedom of the polymer chain 
are at equilibrium. 



polymer chain by a single spring (thereby obtaining a 
Markovian problem) [10] [26] [27], or by assuming that 
the distribution of position of the non-reactive monomers 
instantaneously reaches a local equilibrium assumption 
[8j [9] . These theories were formulated in the context of 
intramolecular reactions, but can be generalized to the 
case of intermolecular reactions as well. Other theoretical 
approaches include the use of the renormalization group 
theory [TTJ[T2] that provides for infinitely long chains per- 
turbative results in the parameter e = 4 — where d is 
the space dimension. Because Markovian theories have 
been used in the analysis of recent experimental works 
on hairpin formation or the folding of polypeptide chains 
[3 [SHI] i and because Markovian approximations (such as 
the quasi independent intervals approximation [28 ) play 
an important role in the study of general stochastic pro- 
cesses, it is important to establish the validity regime of 
Markovian theories and the order of magnitude of non- 
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Markovian effects. 

In a recent work, we proposed another approach of the 
problem, in which the non- Markovian effects are explic- 
itly taken into account by determining the statistics of 
the polymer conformations at the very instant of the re- 
action [29] . This non- Markovian theory predicts that the 
polymer is elongated on average at the instant of reac- 
tion, as can be seen in Fig. |2|a). This elongation does not 
exist in equilibrium conformations [Fig. |2jb)], which are 
assumed to be the reactive conformations in a Markovian 
approach. Due to the elongation of the reactive confor- 
mations, in the non-Mar kovian description the polymer 
centre of mass therefore needs to approach the target less 
closely than it does in the Markovian theory, which leads 
to faster reaction kinetics. 

The main goal of the present paper is to complete the 
initial presentation of this non-Markovian theory of poly- 
mer reaction kinetics and to present new results in the 
case of intermolecular reactions. In particular, we use the 
theory to estimate the magnitude of the non-Markovian 
effects. We find that, for a polymer in three dimensions 
(3D), the non-Markovian effects on the reaction time are 
of the same order of magnitude as the expression of the 
reaction time obtained in the Markovian theory. One of 
the most striking results of the present study is that in 
a one-dimensional (ID) space the non-Markovian effects 
are much stronger: the physics of the diffusion controlled 
reaction in ID is not properly described by a Markovian 
theory, and for long chains the reaction time predicted 
by the non-Markovian theory can be orders of magni- 
tude smaller than in the Markovian approximation. In 
this paper, we also provide an analytical description of 
the average reactive conformation of the polymer, and we 
describe the non-Markovian theory in detail. We com- 
plete the study by showing that it is possible to derive 
systematically scaling expressions for the reaction time 
that take into account the behavior of the monomer mo- 
tion at various time scales. The present paper deals with 
intermolecular reactions and will be completed by an- 
other paper focused on intermolecular reactions such as 
cyclization [30] . 

The outline of this paper is as follows. In the section 
PH we briefly introduce the Rouse model of a polymer 
chain, and we define the notations that we use. Then, 
in the section |III[ we introduce a systematic manner to 
derive scaling expressions for the reaction time, in var- 
ious regimes both in one dimensional and three dimen- 
sional spaces. Afterwards, we present a detailed descrip- 
tion of the non-Markovian theory in the case of a chain 
evolv ing in a one-dimensional space (sections |IV A and 
IV B). In this theory, we explain how to determine the 



a. non-equilibrium reactive conformations 



statistics of the reactive conformations by using a Gaus- 
sian approximation. Writing the equations requires the 
derivation of projection formulas (that describe the av- 
erage and variance of a monomer position given that the 
reactive monomer is at a fixed position) and propagation 
formulas (that describe how the average and variance of 
the monomer position evolve with time). We carry out 




b. Equilibrium conformations in the reactive region 






FIG. 2. (color online) Examples of conformations of a Rouse 
polymer at the instant where the first monomer hits for the 
first time a reactive sphere. We show both non-equilibrium 
conformations (a) that are predicted by the non-Markovian 
theory, and equilibrium conformations (b). The color scale 
codes for the position of the monomers in the chain, and 
the first monomer is assumed to hit the reactive sphere at 
the north (the vertical direction is indicated by a thin black 
line). As shown in the paper, the fact that at the instant of 
reaction the polymer is more elongated than for an equilib- 
rium conformation implies a faster reaction kinetics. In 3D, 
non-Markovian effects on the mean reaction time quantita- 
tive. Importantly, in ID, the non-Markovian effects are much 
stronger and change the scaling relations for the reaction time. 
Parameters: N = 300 monomers and the size of the reactive 
region is a — O.lV^iV/o, with lo the size of a single bond (see 
text). 



a precise comparison of the non-Markovian theory with 

Then, we show 



stochastic simulations in section IV C 



that the non-Markovian and the Markovian theories pre- 
dict different scaling relations for the reaction time as 
a function of the initial distance between the reactants 



(section IV D), and we give analytical expressions that 
characterize the reactive shape of the polymer in section 



IV E Afterwards, we show how to adapt the formulas to 
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the case of a three-dimensional space (section [V}, where 
we also compare the theoretical predictions with simu- 
lations and derive analytical formulas that describe the 
reactive shape of the polymer in various limiting cases. 
We complete the study by considering the effect of the 
position of the reactive monomers in the chain on the 
reaction kinetics (section VD). 



II. THE ROUSE POLYMER CHAIN: 
DEFINITIONS AND NOTATIONS 

We consider the classical model of a Rouse chain of 
N monomers connected by linear springs of stiffness k. 
The monomers experience a frictional drag of coefficient 
£ and diffuse with a diffusion coefficient D = fc#Tb/C 
in the force- field created by their neighbors, with To the 
temperature. Even if this minimal model neglects both 
hydrodynamic interactions and excluded volume effects, 
it captures some of the main features of polymer dynam- 
ics [T7J [18] . Its simplicity makes it suitable to examine 
precisely the different theories of polymer reaction kinet- 
ics, that are in fact non trivial [25]. We denote the micro- 
scopic time scale by To = (/k : which is the typical relax- 
ation time of a bond in the polymer, and the microscopic 
length by Iq = ^/kBT/k^ which is the typical length of 
a bond. We introduce the positions r^, i G {1,...,7V} 
of the N monomers, where quantities in bold stand for 
vectors in the <i-dimensional space. The evolution of the 
probability P(i*i, rjv? t) to find the polymer chain in a 
given configuration at time t satisfies the Fokker-Planck 
equation [T71 \TE\ I3T] : 



dP 
~dt 



N N 



(1) 



where = d/dvi is the nabla operator for the position 
of the i th monomer, and is the force acting on the i th 
monomer. As the monomers are connected by springs, 
this force is related to the monomer positions by a linear 
relation: 



N 



(2) 



where the connectivity matrix M reads: 

/ 1 -1 .. .. 
-12-10 .. 
0-12-1 .. 



M 



\ 



0-12-1 
.. 0-11/ 



(3) 



It is useful to consider the eigenvalues and eigenvectors 
of M because it will enable the definition of the Rouse 
modes, which considerably simplify the description of 
the dynamics of the polymer. Because M is tridiago- 
nal positive symmetric, it can be diagonalized: we write 



M = Q diag(Ai, Aat) where diag(Ai, Ajv) is 

the diagonal matrix with the eigenvalues A^ on the diag- 
onal, and Q is an orthogonal matrix that is normalized 
such that its inverse is its transpose: Q~ l = Q l . The pos- 
itive eigenvalues A^ and the coefficients Qij of the transfer 
matrix can be written explicitly: 



A i = 2{l-cos[(i-l)7r/7V]} 



(4) 



^'co S [(i-l/2)(j-l)7r/iV] (5) 



N 

where 1 < z, j < N , and 5^ is the Kronecker delta sym- 
bol. The definition of the transfer matrix enables us 
to define the Rouse modes a^ with the two (equivalent) 
equations: 



N 



N 



1 Qj 



(6) 



The evolution of the probability P(ai, ajy, t) of ob- 
serving the Rouse modes at time t is given by a Fokker- 
Planck equation that can be deduced from Eq. 0: 



dP 



TO 



This equation shows that, if the modes a^ are inde- 
pendent at some time to, they remain independent at 
all later times t > to. Note that the first eigenvalue 
vanishes (Ai = 0) ; the corresponding eigenmode ai is 
therefore diffusive, it is indeed proportional to the posi- 
tion of the polymer center-of-mass r cm , which is given 
by: r cm = ai/o/V^V- The smallest non-zero eigen- 
value A2 is inversely proportional to the largest relax- 
ation time of the internal conformations of the chain. 
This time is named the Rouse time and is given by 
t r = ^0/^2 — N 2 ro/7r 2 . This mode describes the dy- 
namics of the chain at the length scale loVN. From Eq. 
Q, we also note that the largest eigenvalue is A at — 4, 
the smallest time scale of the internal degrees of freedom 
of the chain is simply To/ '4: it remains of the order of the 
individual bond relaxation time, and the larger modes 
describe the dynamics of the chain at the microscopic 
length scale Iq. For infinite TV, the decomposition of the 
positions into Rouse modes a.j is equivalent to taking 
the Fourier transform of r(s), where s is the curvilinear 
coordinate of a monomer in the chain [T8] . 

In this paper, we are interested in intermolecular reac- 
tions, and we focus on the reaction between a monomer 
(located at position of index p in the chain) and a fixed 
external target of size a. The reaction is assumed to take 
place in a large confining volume V. The p th monomer 
is called the reactive monomer, and we note R = r p its 
position. The position of the reactive monomer R can 
be expressed as a linear sum of the Rouse modes of the 
chain: 



N 



R 



<6|a> 



(8) 
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where the coefficients b{ are easily identified by consider- 
ing Eq. (pj): 



(9) 



In Eq. (J8|, we also introduced the notation |a) that rep- 
resents the full polymer conformation (ai, 8ln)- The 
notation \u) represents the A/'-components column vec- 
tor (ui, ujsf). The quantity (u\ is the transpose vec- 
tor of \u), and for any symmetric matrix A, we define 
= - =1 UiAijVj. Note that quantities in bold 
represent vector in the physical d-dimensional space, to 
be distinguished from the N components vectors noted 
\u). 

We assume that the reaction between the reactive 
monomer and the external fixed reactant is fully trans- 
port controlled and takes place instantaneously as soon 
as the two reactants become closer than a certain cap- 
ture radius a. Then, the reaction kinetics is quantified 
by the mean time T for the reactive monomer to reach 
a sphere of radius a around the external fixed reactive 
site. This time depends on initial conditions chosen for 
the polymer. We chose to study the case where the poly- 
mer is initially at equilibrium, with the condition that 
the initial position of the reactive monomer is Ro- T is 
also called the reaction time. The reaction takes place in 
a confining volume V that is assumed to be large. In par- 
ticular, its diameter V x l d is assumed to be much larger 
than the polymer size (/oV^)- All the theory presented 
in this paper aims at giving an estimate of the reaction 
time T that takes into account non-Markovian effects. 



III. SCALING RELATIONS FOR THE 
REACTION TIME 

A. The root-mean square displacement at different 
time scales 

Before we present the full formalism of the non- 
Markovian theory, we give some qualitative arguments 
that enable the determination of scaling relations for the 
reaction time. We define the important function ip(t), 
that characterizes the stochastic process R(t) in the ab- 
sence of confinement. Assume that at t = 0, the reac- 
tive monomer position is known to be Ro, and that all 
the internal degrees of freedom of the chain are at equi- 
librium. The mean square displacement of R(t) at later 
times t > is called ip(t)/d (where d is the spatial dimen- 
sion). Hence, ip(t) is defined such that as the variance of 
any of the coordinates (X, Y or Z in 3D) of the reactive 
monomer position, given that initially the polymer is at 
equilibrium, and that the initial position of the reactive 
monomer is known: 



ip(t) = var(X(t)|{stat,X },t 

= i N r + 2 & 2 ( 1 - e " 

j>2 



= 0) 

Xjt / r °)/\j. 



(10) 



Here, var(A|£?) represents the variance of the variable A 
given the event B. The expression (10) will be justified 
later in the paper [see Eq. ([51])]. Note that tjj(t) does not 
depend on the particular initial position Ro, but that 
it is a function that involves many time scales A" 1 that 
come from the presence of all the Rouse modes. The 
contribution of each mode to ip(t) is proportional to b 2 , 
and we remind that, by Eq. (pi), b{ can be seen as the 
projection of R on the mode a^. 

From Eq. (10), we can extract the behavior of ip(t) 
at long and short time scales (and at the corresponding 
length scales): 



2Dt 



ift<r (A J R<Z ) 



2{D/N)t if i>r iV 2 (Ai?>WA0 



, (11) 



where we noted AR ~ ^(t) 1 / 2 the typical excursion dis- 
tance by which a monomer moves up to time t. At very 
short time scales, the reactive monomer diffuses as if it 
were disconnected from the rest of the polymer. This 
regime holds as long as the motion occurs at length scales 
AR that remain smaller than the bond length Iq. At 
long time scales, there is another diffusive regime, the 
reactive monomer diffuses with the same diffusion coef- 
ficient as the polymer center-of-mass D cm = D/N. This 
regime holds when the typical distances are large com- 
pared to the polymer size Iq\/N : or, equivalently, at time 
scales larger than the largest internal relaxation time 
tr ~ N 2 tq. At intermediate time scales, all the inter- 
nal time scales contribute to the motion, and it is known 
that the monomer motion becomes sub-diffusive [17] (see 
also Appendix [A]): 



m ~ k I 2 (t/r ) 1/2 



if r <C t < N 2 r 

(Z < AR <C I q Vn) (12) 



Here, n is a numerical coefficient that depends on the 
position of the monomer in the chain: n = for a 

monomer located at the end of a polymer, and n = 
for a monomer located in the interior (see Appendix [A] 
and other references [17]). The smaller value of n for an 
interior monomer is due to the fact that in this case, the 
motion is slowed down by two branches of polymer that 
are surrounding the reactive monomer, instead of only 
one branch for an exterior monomer. Equation (12) in- 



dicates that the motion is subdiffusive, and enables us to 
define an effective walk dimension [32 d w with the rela- 
tion ip ~ t 2 / d ™^ leading to d w = 4 at these intermediate 
length and time scales. 



B. Scaling laws for the reaction time in 3D 

Let us now use the different expressions of ip(t) at the 
various time and length scales in order to derive scal- 
ing relations for the reaction time T in 3D. Given that 
R is diffusive at long times, we expect that when the 
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initial distance between the reactants increases, the re- 
action time reaches a saturating value which is equal to 
the reaction time average over all initial positions in the 
confining volume. We therefore assume that Ro — » oo 
and we discuss the different regimes with the value of 
the capture radius a. In the following, we use two re- 
sults, that are exact for Markovian variables. First, the 
time for a diffusive walker (with diffusion coefficient D) 
in 3D to reach a target of size a in a confining volume V 
is given by the formula: T = V/(47rDa) [33-35 . Second, 
the time needed for a walker that has a walk dimension 
d w > d to reach a target of size a starting from an initial 
distance ro is approximately given by: T ~ (ro — a) dw ~ d V 
[36j [37] . Let us first discuss the case of a large capture 
radius (a ^> Iq^/N). In this case, only the large length 
scales are involved in the process of finding the reactive 
region. At these length scales, by Eq. (11), R(t) behaves 



as a diffusive walker (with the diffusion coefficient equal 
to that of the center-of-mass D/N), and the reaction time 
is therefore given by: 



V 



4n(D/N)a 



(if a>/ ViV) 



(13) 



Let us now assume that the size of the reactive region 
lies in the intermediate regime: Iq <C a <C IqVN. Then, 
the reaction occurs in two steps. The first step con- 
sists in reaching for the first time a region of size loVX 
around the reactive zone. This step is done by diffu- 
sion, with diffusion coefficient D/N, and lasts a time 
T ~ V/(4tt(D/N)Io\^N). The second step consists in 
reaching the reactive region of size a, while the initial 
distance between the reactants is loV~N. At these length 
scales, the motion is subdiffusive [Eq. ( 12 )] and therefore 
this step lasts a time T ~ (loVN -a) dw ~ d V, with d w = 4 
and d = 3. Hence, all together, the reaction time in this 
regime reads: 



V 



V(l VN-a) 



47t(D/N)1 VN 



Dll 



(if/ <Ca<C/ V^), (14) 



where we have added the appropriate microscopic length 
and time scales to obtain an homogeneous formula. 

The last case is that of a very small reactive region 
(a <C lo)- The first step of the reaction still consists 
in reaching the radius lo\/N by diffusion (with diffusion 
coefficient D/N). The second consists in reaching the 
radius Iq by subdiffusion, starting from an initial distance 
loy/N, and the last step consists in reaching the size a 
by diffusion (with the same diffusion coefficient D as a 
single monomer). Hence, the reaction time is a sum of 
three terms: 



V 



V(loVN-l ) 



V 



4tt(D/N)1 VN 



Dll 



4nDa 

(ifa<Z ) (15) 



Finally, we can simplify the equations ( 13 ),( 14),( 15 ) for 
N > 1 to obtain:: 



'VN/(4irDa) 
VVNt /1 3 



V 



Nr /l% + l/{4nDa) 



if loVN <C a 

if l < a < loVN 

if a <C lo 

(16) 



From the last line of equation (16), we get the inter- 



esting observation that in the regime a <C Zo? the reac- 
tion time is the sum of two different times. The part 
T ~ V/(Da), that comes from the diffusive behavior of 
the motion at short time/length scales, dominates only 
for very small length scales (a <C lo/y/N). In fact, in 
the regime lo/VN <CflC ZoV^V, the reaction time does 
not depend on the capture radius a. This fact is a conse- 
quence of the subdiffusive behavior of the motion at the 
intermediate length scales, which implies that the spa- 
tial exploration is compact. It was already known from 
the early analyses of De Gennes [13] and Doi [27] (in the 
context of cyclization) , or with the renormalization group 
theory [TTJ[T2]. However, we are not aware of any existing 
systematic method to derive systematically all the inter- 
mediate scaling laws (16) that appear for intermediate 



values of the reactive sizes. 



C. Scaling laws for the reaction time in ID 

Let us now derive scaling expressions for the reaction 
time in the case of a one-dimensional space. In this case, 
the size of the reactive region can be taken to be a = 
0, the reactive region is simply the point at the origin 
of the spatial coordinate X. The fact that the reaction 
takes place in a volume V means that there is a reflecting 
wall at the coordinate L (in fact, the effective volume is 
V = 2L). The initial distance between the reactants 
is Xq, and the regimes of reactions are determined by 
discussing with the value of Xq. 

Let us first consider the case Xq <C lo- At these small 
length scales, the reactive monomer behaves as if it were 
alone ; it diffuses with a diffusion coefficient D, and there- 
fore the reaction time is: 



LXp 
D 



(ifX <Z ). 



(17) 



The second case is that of a larger initial distance (lo <C 
X <IC loVN). Then, the reactive monomer needs to 
reach the size m a subdiffusive way, and then diffuses 
until it reaches the target by diffusion. The reaction time 
is therefore a sum of two contributions: 



Llo 
D 



l(x - i y 

Dll 



LXl 
Dll 



(if Jo <.X <:l VN). 

(18) 



The last case to consider is Xq ^> lovN, in which case 
the reaction first consists in reaching a distance IqVN 
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from the target (by diffusing with the polymer center-of- 
mass diffusion coefficient D/N), followed by a subdiffu- 
sive step to reach the length an d a diffusive step to 
reach the reactive region: 



D 



L(l VN - l ) 3 L(X - kVN) 



Dll 



D/N 
(if X ^>1 Vn). 



(19) 



At this stage, we have proposed a simple way to de- 
rive scaling arguments for the reaction time both in ID 
and 3D by taking into account the various limiting be- 
havior of the mean square displacement function ip(t) at 
different time scales. Further analysis is necessary for 
various reasons. First, all the numerical coefficients that 
appear in the scaling relations are unknown, which does 
not facilitate the comparison with numerical simulations. 
Second, the derivation of these scaling laws is based on 
scaling arguments that are valid for scale invariant pro- 
cesses, which is not the case here. In particular, the de- 
composition of the reaction between different substeps is 
not obvious. For example, in the 3D case, if the monomer 
reaches a distance loVN from the target, it has a prob- 
ability to escape at distances much larger loy/N before 
it reaches the target, and therefore one could guess that 
subsequent steps of the reaction involve the behavior of 
ip at length scales larger that IqVN. Last, this simple 
analysis is based on arguments that are valid for Marko- 
vian processes, whereas the stochastic process R(t) is 
non-Markovian. Actually, we will find that a more re- 
fined Markovian theory, based on a Wilemski-Fixman 
type approximation, predicts a scaling relation different 
from ( 18 ). There is therefore an ambiguity on what is the 
expression of the reaction time predicted with Markovian 
assumptions in the regime (18). As we shall see below, 



the non-Markovian prediction for the mean reaction time 
is equal to the scaling (18) and is therefore different from 



the Markovian prediction. 

In the next sections, we describe in detail a non- 
Markovian theory that enable a precise determination of 
the reaction time. In the rest of the paper, we choose the 
microscopic length of a bond Iq as the unit of length, the 
typical relaxation time tq is chosen as the unit of time, 
and the unit of energy is ksT. Therefore, we can write 
tq = 1, Iq = 1 and fc#T = 1 in the theory. 



IV. THE NON-MARKOVIAN THEORY IN ID 

A. The renewal equation and the distribution of 
reactive conformations tv 

After these scaling arguments, we present a complete 
theory that enables the precise determination of the reac- 
tion time in this non-Markovian problem. For simplicity, 
we present the complete theory in ID. In the generaliza- 
tion of the theory to a 3-dimensional space, geometrical 
effects appear and will be described in section (V) We 



stress that, even if the ID case is quite artificial in the 
context of polymers, it could be relevant in other contexts 
such as the study of first passage properties of a noisy 
moving interface [38 j. In ID, the observable R is iden- 
tified with its first coordinate X, and we are looking for 
the time for the reactive monomer to reach the position 
Xf = 0, while its initial position was Xq > with a equi- 
librium configuration for the rest of the chain. While the 
dynamics of the position of the reactive monomer X(t) 
is non Markovian, the evolution of the full polymer con- 
formation | a) is Markovian and obeys a renewal equation 
[3T] which is the starting point of our analysis. Let us 
consider a polymer conformation \a) such that (b\a) = Xf 
(i.e., such that the reactive monomer is inside the reac- 
tive region at position Xf). Consider the situation where 
the polymer does not react when X reaches the position 
Xf. Then, observing such a conformation \a) at time 
t necessarily implies that the polymer has reached the 
target for the first time at a time t' < t, with some con- 
formation \a r ). Therefore, if we define f(\a'),t') as the 
probability density that, starting from the initial distri- 
bution, the reactive region is reached for the first time 
at t' with a configuration la 7 ), we can write the following 
renewal equation: 

P(|a),t|{stat,X },0) = 



dt' / d\a!) f(\a'),t'\)P(\a),t-t'\ |a'),0). (20) 



Here, d\a) = dai...dajsf, and P(\a), t|{stat, X }, 0) is the 
probability of observing a configuration \a) at t in the 
absence of reaction when the initial distribution at t = 
is an equilibrium distribution with the reactive monomer 
in position Xq. Similarly, P(\a),t — 1'\ |a'),0) is the prob- 
ability of observing the configuration \a) at t given that 
the configuration \a') was observed at t = 0. We in- 
troduce the splitting probability distribution 7r(|a)) = 
7r(ai, cln) that represents the probability density of 
observing a configuration \a) when the reaction takes 
place: 



7r(\a)) 



f 

Jo 



dt f(\a),t). 



(21) 



This splitting probability distribution depends on the ini- 
tial conditions, but we do not make it appear explic- 
itly in the notations for simplicity. Taking the Laplace 
transform of the renewal equation ( 20 ) and expanding for 



small values of the Laplace variable yields the two follow- 
ing relations, valid for all the conformations \a) such that 
(b\a) = X f . 



J d\a) tt(|o» = 1, 
T(X )P stat (\a)) = 



(22) 



dt 



P(\a),t\ir,0) - P(\a),t\{st&t,X o },0) . (23) 
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Here, we have introduced P s tat(|&))> that represents the 
probability of observing a given configuration in the sta- 
tionary state. The quantity 1 1 tt, 0) is the probabil- 
ity of a configuration \a) at t given that the configuration 
at t = is taken from the splitting probability n, and it 
is given by: 

P(|o>,t|7r,0) = / d\a')n(\a'))P(\a),t \ \a'),0). (24) 



The equations ( 23|24 ) together with the normalization 
condition (22) form an integral equation that completely 



defines the splitting probability tt and the mean first re- 
action time T, but which is very difficult to solve in the 
general case. From Eq. (23), we can derive several sets 
of equations that must be satisfied and that link the re- 
action time to the moments of the splitting probability. 
First, we need to reinterpret Eq. (23 ) (which is valid only 
for configurations such that (b\a) = Xf). We note that 
the probability density to observe the reactive monomer 
in position Xf given that the rest of the polymer has 
the conformation \a) is simply given by S(Xf — (b\a)). 
Therefore, using the Baye's formula, we can write: 



P(\a))S((b\a)-X f )=P(X f )P(\a)\X f ). 



(25) 



Hence, multiplying the integral equation (23) by 5 (Xj — 
(b\a)) and using the trick (25) enables us to write (23) in 
a slightly different way: 

TP stat (X f )P st&t (\a)\X f ) = 



dt 



P(Xf,t\ir,0)P(\a),t\X f ,t;ir,0) 



- P(X f , t|{stat, X }, 0)P(\a),t\X f , t; {stat, X }, 0) 

(26) 

where P(\a),t\Xf,t]ir,t) is the probability of observing 
the configuration \a) at t given that X = Xf at the same 
time t and that the distribution of modes at t = was 
tt. Note that Xf = 0, but we keep the notation Xf so 
that there is no confusion with the initial time t = 0. 
Similarly, P stat (\a)\Xf) is the stationary probability to 
observe a configuration given that the value of the ob- 
servable is Xf (in the absence of reaction). Now, the 



equation (26) is valid for any value of \a) (not only for 
those that satisfy (b\a) = 0), and it is exact if all the 
quantities are evaluated by taking into account the con- 
fining volume V. At this stage, we do a large volume 
approximation: we assume that all the terms appearing 
in (26) can be replaced by their expression in unbounded 
space, except for the term P stat (Xf), which is replaced 
by the inverse of the confining volume 1/V: 



t(X f )~l/V (V^oo). 



(27) 



We will derive below the large volume asymptotics of 
the mean first-passage time. Note that, in ID, if L is 
the distance that separates the target from the reflecting 
wall, the value of the confining volume is 2L. Noting 



that the distribution P(\a), t\Xf, t; 7r, 0) is normalized to 
1, it is clear that the integration of Eq. (26) over all the 
conformations la) leads to: 



TV~ l = 



dt 



P(X f ,t\ir,0) -P(X f , t\{sta.t,X },0) . (28) 



This expression generalizes the results obtained for 
Markovian systems [36j [33 [39j [40] , and makes it clear 
that the mean first passage time can be expressed as time 
integrals of propagators. The equation (28) is essential 



and is at the basis of all our estimates of the reaction 
time in this paper. As P(Xf,t\ir,0) is the probability 
of observing the reactive monomer at a position Xf at t 
given that the initial conformational statistics is the split- 
ting probability n, determining the splitting probability 
distribution tt is a key step in determining the kinetics 
of polymer reactions. This step is however highly non- 
trivial and consists (in principle) in solving the integral 
equation (23). Because this integral equation involves 



functions of N variables, its solution seems out of reach 
with analytical tools. A natural attempt to overcome 
this difficulty is to assume that the splitting probabil- 
ity tt can be replaced by the stationary probability of 
conformations restricted to conformations \a) such that 
(b\a)=Xf. 

tt(\o)) ~ P s tat(|a)|X = Xf) (Markovian approx.) (29) 

We call this approximation the Markovian approxima- 
tion: all the memory effects are neglected since it is as- 
sumed that the polymer reaches instantaneously its equi- 
librium distribution. In particular, the dependance of the 
splitting distribution with the initial conditions cannot 
be addressed in this approximation. As shown elsewhere 
[29} [30] , the corresponding approximation in the case of 
intramolecular reactions gives the same results as the 
classical Wilemski-Fixman approximation with a certain 
choice of sink-function [SJ [9j [25] . Using the Markovian 



approximation ( 29 ) leads to the following approximation 



of the first propagator appearing in Equation (28): 



P(X f ,t\n,0) 



1 



1/2 



exp 



(Xf? 



(30) 



where the function ip is given in Eq. ( |10| ). This leads to 
the Markovian estimate of the reaction time: 



^Markovian V 



f 

JO 



dt 



V2tY$ _ 



1 — exp 



2ip 



■ (31) 



However, as shown below, this Markovian approximation 
does not support the comparison with numerical simula- 
tions. We show now a method to go beyond this Marko- 
vian estimation of the reaction time. 



The general equation (26) does not only lead to the 



estimate of the reaction time: other relations can be de- 



rived. For example, multiplying Eq. (26) by ai and in- 



tegrating over all the modes leads to the definition of 
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TABLE I. Summary of the notations used in the non-Markovian theory in ID, with references to the equations where the 
quantities are defined. 



ml 

im 
nr (t) 

mf at 

{stat.Xn } 

4 stat ' Xo} (t) 

{stat,X },X f 



average of a% at the reaction 
average position of the i th monomer at the reaction 



average of a% at a time t after the reaction [Eq. (37)] 



average of a% at a time t after the reaction given that X = Xf at t [Eq. (42)] 
average of X at a time t after the reaction [Eq. (39)] 
covariance of a% , dj at the reaction 



covariance of a*, dj at a time t after the reaction [Eq. (38)] 



covariance of a», aj at a time t after the reaction given that X = at t [Eq. (43)] 
covariance of X at a time t after the reaction [Eq. (40)] 
average of a% at equilibrium (for i > 2) 



stat,* / stat.Xn 
{stat,Xo},Xj , x 

7ij l r J 



m 

Si i/N 



average of ai given that X — Xo and that the polymer is at equilibrium [Eq. ( |47| 
average of a* at £ given that, at £ = 0, one has X — Xo and the polymer is at equilibrium [Eq. (49)] 
average of a% at £ given that X = Xf at £, and that at £ = the polymer is at equilibrium 

with a reactive monomer in position X — Xo [Eq. (53)] 
covariance of a*, dj at equilibrium (for i, j > 2) 



) covariance of a^, dj given that X — Xq and that the polymer is at equilibrium [Eq. (48)] 



covariance of dj at £ given that, at £ = 0, one has X = Xo and the polymer is at equilibrium [Eq. (50)] 
covariance of dj at t given that X — Xf at t, and that at t = the polymer is at equilibrium 



with a reactive monomer at position X — Xq [Eq. (54)] [it is also noted 7, 



{statjXo},* 



(*)] 



covariance of X at t given that, at t = 0, one has X = Xo and the polymer is at equilibrium [Eqs. (51 10)] 
position of the 2 th monomer in the chain (0 < Si < 1) 



another set of necessary conditions on tt: 



dt 



'p(X f ,t\ir,0)fi' X ' 

P(X / ,t|{stat,Xo},X / )/i 2 {stat ' Xo} '^l =0, (32) 



where /i^' f is the mean value of ai at t given that X = 
Xf at t and that the initial distribution at t = is the 
splitting distribution tt. Similarly, ^{ stat ' X °}' X / j s ^ ne 
mean value of at time t given that X = Xf at t and 
that the polymer was in a stationary state at t = with 
an initial reactive monomer position X = Xq. Note that 
all the notations of this section have been summarized in 
the table H 

Another set of equations is obtained in a similar way: 
multiplying Eq. (26) by aja^ integrating it over all the 



modes and using Eq. (28) leads to: 



C°° r / \ 

/ dt \P{X,,t\v, 0) ( 7 J* + ft ' ' - "T ) 



P(X /5 £|{stat,X },0)x 

{stat,X },* , {stat,Xo},Xj {stat,Xo},X/ 



stat,* 



= 0, 
(33) 



where 7^'* is the covariance between and aj at t given 
that X = Xf at t and that the initial distribution at 
t = is the splitting distribution 7r, while ^ stat ' x °}'* ^ s 



the covariance of a^, aj at t given that X — Xf at £ and 
that the initial distribution at t = is the stationary dis- 
tribution with the reactive monomer located at Xq. The 
quantity cr^ at '* is the covariance of a^a^ at the station- 
ary state, given that X = Xq. The three sets of equations 
( 28|32|33 ) must be satisfied and they provide constraints 
on the possible forms of the splitting distribution tt. 

We now make the key-hypothesis of the non-Markovian 
theory: we assume that the splitting distribution 7r(|a)) 
can be accurately described by a multivariate Gaussian 
distribution. The distribution tt is therefore fully char- 
acterized by the averages (which we call m*) and the 
covariance matrix (denoted afj) of the variables a^. All 
the multivariate Gaussian distributions are not good can- 
didates for the splitting probability tt: the moments 
and (jj, must reflect the fact that X is known with cer- 
tainty to be Xf for this distribution, which implies the 
two conditions: 



(b\m")=X f ; <x» = |0>. 



(34) 



In particular, Eq. (34) states that the matrix a 71 " can- 
not be inverted: this is due to the fact that tt is pro- 
portional to a delta function: 7r(|a)) ~ S((b\a)), and the 
distribution tt is a generalized multivariate Gaussian as 
discussed for example by Eaton [41]. Besides satisfying 
the conditions (34), th e momen ts of tt must satisfy the 
three sets of equations ( 28|32|33 ), which are then used as 
a closed system of self-consistent equations that allow to 
calculate the moments m 71 ", afj and the reaction time T. 
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The last step needed to characterize the non Markovian 
theory therefore consists in relating the moments of the 
splitting distribution raf and af 7 - to the quantities /J^ ,Xf 

32|33|. 



and 7^'* that appear in Eqs. 



(28 



This will be 



done through propagation and projection formulas, that 
we derive now. 

Let us first describe propagation formulas. We call 
/if (£) and jfj (t) the average and covariance of the modes 
di at t given the splitting distribution at t = 0. It is 
well known that the Fokker-Planck equation ^ admits 
Gaussian solutions and that the evolution of /if and 7^ 
satisfies the following equations [31]: 



dt 



(35) 

(X i + X j )^+25 ij . (36) 



The actual values of /if (t) and 7^-(t) are found by solv- 
ing |35|,(36) with the initial conditions /if (0) = mf and 

75(0) 



erf,. We find: 



m i e 

% (1 



-2\it\ 



/A, 



(37) 
(38) 



These formulas describe how the mean vector and the 
covariance matrix are modified with time, and we call 



them "propagation formulas". Note that Eq. (38) is 
written with the convention that (1 — e~ 2Xlt )/Xi = 2t 
(withAi=0). 

We define X n (t) and Vv(^) as the average (and vari- 
ance) of X at given that the initial distribution of the 
monomers was 7r. Because X = (b\a), these two quanti- 
ties are simply given by: 



eE(X,£|tt,0) = 
var(X,£|7r,0) 



= (b\^\b). 



(39) 
(40) 



Then, we can write the expression of the distribution of 
X at t given ir at t = 0: 



P(X,t\n,0) 



1 



1/2 



exp 



2 Mt) 



(41) 



We now derive projection formulas. An explicit expres- 
sion for /J^ ,Xf (the mean of ai at t given a particular value 
of Xf at the same time t) can be found by adapting the 
formulas on conditional probabilities given for example 
by Eaton [41J (Appendix p} : 



Mi 



17*16) 
(b\r\b) 



((b\f)-X f ). (42) 



Here, |e$) represents the i th basis vector (all its elements 
are except for the i th which takes the value 1). Equation 
( |42| ) is in fact very general and states that, if N variables 
di have a mean vector /if and a covariance matrix 7^, 
then the average of over all configurations \a) such that 
(b I a) = Xf is given by the relation (42). 



A similar calculation leads to a second projection for- 
mula for the covariance matrix, that enables to calculate 
the covariance of a^, aj at t given that the position Xf is 
known at the same time t: 



1, 



ir,X f 



(e t \r\b)(e 3 \r\b) 
(b\r\b) 



7f 



(43) 



Note that j njX does not depend on the value of Xf, 
which is why we just note it 7^'* instead of 7 7r ' x z, to the 
difference of \±*' Xf which depends linearly on the value 
of Xf. We call the equations (42) and (43) "projection 
formulas": they describe how the mean and the covari- 
ance of the modes are modified when one restricts the 
modes to be on the hyperplane of equation (b\a) = Xf. 
We now describe how choosing the initial moments 
stat,x an( j a s t at >* ^hat j-^e initial value of the ob- 
servable is Xq, the other degrees of freedom being at 
stationary state. Let us temporarily assume that Ai > 0. 
Then, at stationary state, the moments of a$ are equal 
to: 



; al 



stat 



Applying the projection formulas ( [42} 43), we get: 
X h 



stat,^n 

m- 



stat,* stat 
ij = a ij 



Si 



Ai(6|o- stat |6)' 

(e i |o- stat |6)(e i |o- stat |6) 



(b\a stat \b) 



AiAj<6|a 8tat |6)' 



(44) 



(45) 



(46) 



These formulas are valid under the hypothesis that Ai > 
0. Taking the limit Ai — » leads to: 



stat,J^o 



rn 



stat,* 



SnXo/h 



if i,j > 2 
if 3 > 2, i = 



(47) 
(48) 



Applyi ng the propagation formulas ( 37|38 ) to equations 
fl47p8l ) leads to: 



/i, 



{stat,Xo} 



stat,X e -A 4 t 



{stat,Xo} 



% (1 



2A,t) /A 



(49) 



e _ Ai t e _ A . V stat,*_ 



(50) 



From (49), we deduce that X remains on average at the 



position Xq, while the expression (50) enables us to ex- 
plicitly calculate function iJj: 



N 



t/>(t) = var(X,£|{stat,X },0) = ^ hb^f^^ (51) 



Reporting the equations ( 5Q|48 ) in this definition leads 
to the expression of ip(t) that we had given earlier [Eq. 
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(10)]. We now derive the projected quantities. From 
equations ( 37|38 ), we get, for i > 2: 



( e .| 7 {stat,X }| 6) = 6 . /A . (1 _ e -A,t ) . 



(52) 



Therefore, using the projection formulas (42, 43), we ob- 
tain, for i,j > 2: 



{ s tat,x },x / _ h(l-e Xit ) 



(X ~ X f ) , 



{stat,Xo},* 



A,- 



bibj(l - e- Ait )(l - e-^*) 



(53) 



(54) 



These expressions are valid for z, j > 2, but suffice to fully 
determine the moments of 7r, because these moments also 
satisfy the condition ( [34] ). 

Let us also write the explicit expression for the effective 
propagator P(X, t|{stat, X }, 0): 



P(X,£|{stat,X },0) 



(27H/0 



^ 6XP {" 



(x - x ) 

2V(*) 



(55) 



Last, the general estimate of the reaction time p8[ ) can 
be written more explicitly in the non-Markovian theory: 



T 
V 



dt 



2tt 



>y 2 



exp 



xl 

2^ 



exp 




(56) 



At this stage, we have completely defined the non- 
Markovian theory. This theory consists in determin- 
ing the mean vector \m K ) and covariance matrix a 71 of 
the splitting distribution tt by solving the set of self- 
consistency equations ( 32|33 ) . All the quantities appear- 
ing in these equations can be explicitly related to and 
afj through the propagation formulas (37 38|49|5Q ), the 
projection formu las ( 42|43|53|54 ) and the expressions of 
the propagators ( 41|55 ) . Once the moments of the split- 
ting distribution are determined, one can estimate with 
Eqs. (39, 40) the mean X n (t) and the variance ^(t) of 
the trajectory of the reactive monomer at a time t after 
the first time it reached the reactive position. These two 
quantities appear explicitly in the expression ( 56 ) for the 
reaction time, which can be evaluated. 



B. Simplified non-Markovian theory in ID: the 
stationary covariance approximation 



is a difficult task when N is large. Up to now, the only al- 
ternative to the full non-Markovian theory is the Marko- 
vian approximation (29), in which all memory effects are 



neglected. We propose here another alternative, that we 
call the "stationary covariance approximation" , in which 
the covariance matrix of the splitting distribution is as- 
sumed to be well approximated by the covariance matrix 
of the modes in the stationary state restricted to config- 
urations such that X = X : 

^ r^j —Stat,* /r >7\ 



Note that a^ at '* is given by Eq (48), and that this ap- 



proximation implies that ^(t) ~yj[t). Within this ap- 
proximation, afj cannot be expected to satisfy the set of 
equations (33), which must therefore be released. The 



moments are then determined by solving the remaining 
set of TV — 1 equations (32 ), whose expression can be made 



a little simpler, as we obtain by using (|52[) for i > 2: 



exp 



xl 

2xp 



m- e 



-Xit 



bi(l 



-Xit 



-Xit 



)X 



Xi ip 



Xi ip 



= o, 



(58) 

where X n can be written as a function of the modes 
with i > 2 only: 

N 
i=2 

This expression follows from the fact that (blm^) = 0. 
The equations ( 58|59 ) completely define the moments 
m* . Last, the expression of the reaction time in the 
stationary covariance approximation reads: 

r dt f ( x i\ ( x * 



T 
V 



2i/j 



2ip ) 



(60) 

The stationary covariance approximation can be seen as 
an intermediate theory between the Markovian approx- 
imation and the complete non-Markovian theory: while 
being much simpler than the non-Markovian theory, it 
catches some memory effects by considering the average 
positions of the monomers at the instant of the reaction. 
We remind that, in all these theories, apart from the 
hypotheses on the shape of the splitting probability, we 
made a large volume approximation. In appendix [Cj we 
present the method that we used to obtain the numerical 
integration of the theory. The theoretical results in each 
version of the theory are compared with simulations in 
the next section. 



Before we compare the results of the non-Markovian 
theory with numerical simulations, we propose a simpli- 
fied version of this theory. A simplified version is nec- 
essary because the complete non-Markovian theory re- 
quires to solve a set of N — 1 + (TV — l) 2 equations, which 



C. Comparison with numerical simulations 

We performed stochastic simulations with a simulation 
algorithm that is described in details in the appendix [D) 
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In short, the polymer evolves in a box between X = 
and X = L, the position X = L is reflecting for the 
first monomer, whereas X = is absorbing for the same 
monomer. All the other monomers do not see the bound- 
aries. At each run, one generates an initial equilibrium 
configuration for the polymer, which is then translated 
so that the initial position of the reactive monomer is Xq. 
Then, the polymer evolves at each time step (of fixed size 
At) according to an adaptation of the Brownian dynam- 
ics simulation algorithm of Peters et al. [42] . If it is close 
from the reacting region located at X — 0, at each time 
t, one computes the probability to be absorbed between 
t and t + At, and the generation of a random number 
enables to decide whether or not the simulation has to 
be stopped at this time step. If it is not the case, all the 
monomers evolve under the influence of a Gaussian white 
noise and in the force field of their neighbors, according 
to the Langevin equation that corresponds to the Fokker- 
Planck equations ([!]). One simulation run stops as soon 
as the first monomer reaches the absorbing boundary, 
and one records the position of the other monomers at 
this instant. Although the sampling of the stochastic tra- 
jectories is not exact, we expect it to be more and more 
precise as the time step gets smaller (At — ^ 0). 

On Figure [3J we present the results of numerical sim- 
ulations for N = 20, for which the confining volume is 
about 25 times the size of the polymer (L = 2b\fN) 
and the time step is At = 0.0005, that is 2000 times 
smaller than the relaxation time of a single bond. There 
is a very good agreement between the reaction times pre- 
dicted by the non-Markovian theory and the simulation 
points. The predictions that use the stationary covari- 
ance approximation are slightly less good but the two 
theories differ by only less than 15%. Finally, the predic- 
tions of the Markovian theory are in clear disagreement 
with the simulations, as they differ by a factor of roughly 
2. These remarks remain true for a larger value of TV 
(N = 40), as can be seen in the inset of Fig. |3| where 
the reaction time in the Markovian approximation and 
in the simulations differ by a factor 3, whereas the non- 
Markovian theory predicts correct values of the reaction 
time. For this value of N however, we do not have any 
estimate for the non-Markovian theory, we only have the 
result in the framework of the stationary covariance ap- 
proximation and they are in good agreement with the 
simulations. 

In the numerical simulations, the final positions of the 
monomers are recorded at each run: this enables us to 
determine whether or not the Gaussian approximation, 
which is the key hypothesis of the non-Markovian the- 
ory, is a good approximation. For each value of X , we 
computed the values of the modes at the end of simu- 
lations, rescaled it by the empirical means and variance, 
and plotted the histogram of the obtained distribution. 
Repeating this procedure for all the values of X and all 
the values of z, one obtains several histograms that are 
all superposed in Fig. [4] If the Gaussian approximation 
is good, these histograms should resemble the centered 




FIG. 3. (color online) Reaction time T divided by the confin- 
ing volume V as a function of the initial distance Xq between 
the reactants in ID for N — 20 (main figure) and N — 40 (in- 
set). The reactive monomer is the first monomer. Symbols: 
results of numerical simulations. The half-length of the error- 
bars is equal to twice the standard deviation of the mean of 
the simulated samples, so that the error-bars represent the 
95% confidence intervals. Red line: non-Markovian theory. 
Blue dash-dot line: non-Markovian theory with the station- 
ary covariance approximation. Dashed green line: Markovian 
approximation. The color code are the same for N — 20 
(main figure) and N — 40 (inset). Parameter values: for 
N = 20, there are 39600 simulation runs for each value of Xo, 
with a time step At = 0.0005 in a volume size L = 25V / iV ; 
for N = 40, there are 4000 simulation runs, the time step is 
At = 0.001 and the volume L = The results of the 

complete non-Markovian theory (without the stationary co- 
variance approximation) for N = 40 are difficult to estimate 
and are not shown. 



Gaussian distribution with variance 1, which is also rep- 
resented in Fig. [4j All the marginal distributions of the 
di visually resemble to a Gaussian distribution, both in 
linear scale (Fig. [41 main figure) and semi-logarithmic 
scale (Fig. [4j inset): deviations from Gaussian distri- 
butions are very small. These deviations can be quan- 
tified by the p-values obtained by statistical tests, such 
as the Jarque-Bera test of normality. We found that the 
p- Values obtained for this test are typically of order 1 for 
samples of sizes n ~ 10 4 , whereas the p-values are much 
smaller for larger sample sizes (n c± 10 5 ), meaning that 
deviation from normality cannot be easily detected with 
this statistical test unless the number of samples is larger 
than about 50, 000. 

Then, we compared the theoretical values of and 
gJa to the ones we measured in the simulations. We chose 
to study only the moments of the first mode i = j = 1 
(we remind that a\/y/N is equal to the center-of-mass 
position). On Fig. [5| we represented the values of the 
average position of the center-of-mass position at the in- 
stant of reaction for several values of Xq. For small values 
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a. (centered and rescaled) 



FIG. 4. (color online) Superposition of the empirical probabil- 
ity distributions of all the modes a% (rescaled by their means 
and variance) in ID for N — 20 for all the values of Xo corre- 
sponding to the results shown in Fig. [3] The red curve is the 
normalized centered Gaussian function. Inset: same graph in 
semi- logarithmic scales. This figure shows that the Gaussian 
approximation for the splitting distribution is a good approx- 
imation. 
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FIG. 5. (color online) Average position of the polymer center- 
of-mass at the instant of reaction for N = 20 in ID as a func- 
tion of the initial distance between the reactants Xo. The 
reactive monomer is the first monomer. Red circles: results 
of the numerical simulations (the 95% error-bars are smaller 
than the size of the symbols). Red line: non-Markovian the- 
ory. Dash-dot blue line: stationary covariance approximation. 
Dashed green line: Markovian approximation (the moments 
mj vanish in this approximation). 



of X , there is a good agreement between theory and sim- 
ulations, while for larger values of Xo, the prediction of 
the stationary covariance approximation differs by 9.8% 
from the simulations. When one compares the simula- 



tions with the complete non-Markovian theory, one ob- 
tains a smaller difference of 6.9% for the position of the 
center-of-mass, which is however statistically significant. 
The same remarks hold true for the variance of the po- 
sition of the center-of-mass that is represented in Fig. [6j 
for large values of Xo, there is a difference between theory 
and simulations of about 48% (stationary covariance ap- 
proximation) and 20% (complete non-Markovian theory). 
In the simulations, the position of the first monomer is 
not exactly zero and this causes an uncertainty on the 
empirical value of the center-of-mass position. However, 
in the simulations represented on Figs. |5|6[ the average 
position of the first monomer at the instant of reaction 
is always less than 0.04 and it is not likely that this un- 
certainty can explain the discrepancy between the theory 
and the simulations. Hence, even if the non-Markovian 
theory is much more precise than the Markovian the- 
ory, it does not seem to be an exact theory as it does 
not predict exact values for the moments of the splitting 
probability distribution. We did not expect it to be exact 
anyway. 
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FIG. 6. (color online) Variance of the position of the center-of- 
mass at the instant of reaction for N = 20 in ID corresponding 
to the data of Fig. [3] The reactive monomer is the first 
monomer. Red circles: stochastic simulations. Continous red 
line: non-Markovian theory. Dashed green line: stationary 
covariance approximation (and non-Markovian theory). 

The conclusions to be drawn from this section are the 
following. First, the Gaussian approximation is an excel- 
lent approximation for the splitting distribution. Stan- 
dard normality test such as the Jarque-Bera test cannot 
reject normality of any of the marginal distributions for 
each mode except if the size of the samples is larger than 
50, 000. However, there is no exact agreement between 
the measured values of mf,a? H and the observed ones. 
The theoretical estimate of the reaction time is very pre- 
cise in both the stationary covariance approximation and 
the complete non-Markovian theory. The complete ver- 
sion of the non-Markovian theory and the stationary co- 
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variance approximation give very similar results, and the 
second one can therefore be used in order to obtain very 
precise estimates of the reaction time and the average 
reactive shape of the polymer. 



D. Different scaling relations in the Markovian 
approximation and the non- Markovian theory 



We now focus on the comparison between the differ- 
ences between the predictions of the Markovian and non- 
Mar kovian theories. For simplicity, we restrict the study 
of the non-Markovian theory to the case of the stationary 
covariance approximation, and we assume that the reac- 
tive monomer is the first monomer of the chain. On Fig. 
[7J we have represented the theoretical estimates of the re- 
action time as a function of the initial distance between 
the reactants Xo, both in the stationary covariance ap- 
proximation and the Markovian approximation. In this 
figure, it is clear that both theories predict the same lin- 
ear scaling of T with Xo for both small and large Xo, 
in agreement with the scaling arguments (17) and (19). 



For small Xo, the linear scaling of T with Xo comes from 
the diffusive behavior of the monomer motion at large 
and short time scales, as can be seen from the asymp- 
totics of the function ip(t). As can be observed on Fig. 
[7| the regime of intermediate Xo is remarkable, because 
the Markovian and the non-Markovian theories predict 
very different reaction times in this regime (the predic- 
tions differ by a factor 10 for N = 320 and X = 3.6), 
and the slope of the curves in the log-log plot of Fig. 
[7] are quite different, suggesting that the Markovian and 
non-Markovian theories predict different scaling relations 
in this regime. Furthermore, as can be observed on Fig. 
[8j the maximal ratio of the two estimates of the reaction 
time increases as y/~N. The fact that the difference be- 
tween the two theories can be arbitrarily high for large 
N also suggests that the two theories do not predict the 
same asymptotic relations for T(Xq). The rest of this 
section is devoted to an analytical determination of the 
scaling laws that can appear in both theories. 

The first step of the analysis consists in identifying the 
correct scaling of all the quantities appearing in the equa- 
tions in order to obtain a theory that does not depend on 
N. By Eq. Q, in the limit of large X, the eigenvalues 
are approximated by X q ~ (q — l) 2 7r 2 /X 2 , and from Eq. 

f, we find that b q ~ y / 2/N for q > 2. The fact that 
~ 1/X 2 suggests the definition of the rescaled time 
r = t/N 2 . The theory for infinite X is non-trivial only 
when the parameter Yq = Xq/a/X is fixed as X — >• oo; 
the rescaled initial distance Yq is therefore the initial dis- 
tance between the reactants in the unit defined by the 
typical polymer length y/~N. The correct scaling of the 
moments must leave Eq. (58) invariant with X. We 
find that, if we define M q = m^\/X, Eq. (58) does not 
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FIG. 7. (color online) Rescaled reaction time T/(VN 3/2 ) in 
ID as a function of the rescaled initial distance Yo — Xo / y/N . 
Continuous thick curves: non-Markovian theory (upper red 
curve: X = 80, lower green curve: X = 320). Dashed 
thick curves: Markovian theory (upper red curve: X = 80, 
lower green curve: X = 320). The oblique black line repre- 



sents asymptotic form (65) of the Markovian approximation 
(T = 0.3516 VX 3/2 y 2 )- The oblique dashed blue line repre- 
sents the scaling T ~ Yq, with an arbitrary prefactor. This 
figure shows that the Markovian and non-Markovian theories 
can predict very different values in the regime of intermediate 
initial distances Xo and large X. The reactive monomer is 
the first monomer. 



X 
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FIG. 8. (color online) Maximal ratio T/Tjyiar kovian for several 
values of X, where T is the non-Markovian reaction time in 
the stationary covariance approximation, and Tjyiarkovian the 
Markovian estimate of the reaction time. The oblique dashed 
red line is proportional to X 1//2 . The reactive monomer is the 
first monomer. 
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depend on N anymore, as it reads: 



exp I — I — exp 



2^ 



<2 2 7T 2 ^ 
1 



(61) 



In this equation, Y^{r) is the rescaled reactive trajec- 
tory and ^(t) is the rescaled mean square displacement 
function, which are given by: 



r ) (62) 



9=1 



4(1 - e -9 t r) 



9=1 



q 2 ir 2 



(63) 



Reporting these quantities into the expression for the re- 
action time implies the following asymptotic relation: 



VN 3 / 2 



f 

Jo 



dr 



exp 



Y 2 

_ 7T 

2^ 



exp 



2^ 



(64) 



where the term on the right hand side depends only on 
Yo and not explicitly on N or V. 

For large values of Yb, all the coefficients M q reach 
a fixed asymptotic value and Y n becomes independent 
on Yq. Evaluating (64) by using the large time ap- 



proximation for the integrand leads to the scaling law 
T ~ VN 3 / 2 Yo/2. The asymptotics of the reaction time 
with Yq in this limit is therefore the same for both Marko- 
vian and non-Markovian theories. We now focus on 
the limit of small Yb, where the Markovian and non- 
Markovian theories predict very different values for the 
reaction time. 

The asymptotics of the reaction time T for small Yb 
in the Markovian approximation can be readily found, 
because in this approximation we can write Y n = 0. In- 
serting this equality into Eq. (64) and expanding the 
integrand for small values of Yq leads to: 



T 



Markovian 



VN 3 / 2 



Jo 



dr 



2a/2^ 3 / 2 



0.3516 Y 2 (65) 



integration of the equations for finite N suggest that the 
solutions of the equations have a structure of boundary 
layer when Yq — » 0. The fact that the motion is subdiffu- 
sive at short time scales leads to the definition of the time 
scale to = Y 4 . The function Y n (r) is expected to vary 
at this time scale. The contribution of the first modes 
in Eq. ( [62] ) implies that Y^(r) also varies at the scale 
1. For small Yo, these two time scales are very different, 
which leads us to postulate the following boundary layer 
structure for Y^-(r): 



YAr) 



Y y(T/t ) 
Y [l-ao g(T)} 



T < 1 

r > t 



(66) 



where the parameter ao tends to as Y 0. The rela- 
tion between ao and Yo can be linked to the asymptotic 
form of y(u) and g(r) in the matching region. Let us 
assume the existence of a positive coefficient j3 such that 
g(t) ~ A/t@ for r — >> 0. In this case, the matching condi- 
tion at the intermediate scale to <C r <C 1 imposes that 
y(u — ^ oo ) ~ 1 — A/uP , and that a = t/ Q j = F 4/3 . 

Let us introduce the parameter e that is a matching 
time scale such that to <C e <C 1. By the boundary layer 
hypothesis (66), Y n is well approximated by Yoy(r/to) 
for t < e, whereas it is equal to Yq[1 — aog(r)] for r > 
e. Therefore, we can evaluate the integral appearing in 
Eq. (64) by separating the contributions coming from 
the times r E]0, e[ and r > e. With this procedure we 
obtain the following expression for the reaction time for 
Yq 0: 



VN 3 / 2 



re/to 

Yo 3 

JO 



du 



+ 



Y r 



2+4/3 







2tt 



dr 



9(r) 



[*(r)] 3 / 2 



(67) 



From Eq. (67), it is clear that the scaling of T with 



Yo depends on the coefficient ft: if ft < 1/4, we have 
T ~ Yo 2+4/3 . whereas for ft > 1/4 the scaling is T - Y§ . 
Identifying the coefficient ft is therefore an essential step 
of the theoretical analysis. In the following, we show that 
the only value of ft that is consistent with the theory is 
ft = 1/2. First, we identify the behavior of the moments 
M q when Yq — » that is consistent with the boundary 



layer structure (66). Using Eq. (61), one readily finds 



that the moments M q can be calculated as a function of 



Y n by the formula: 



Note that this integral exists, because for large r we have 

\I/ ~ r, whereas the small r behavior is \I/ ~ r 1//2 (Ap- roc ^ ^_3.^ ^_ q 2 7r 2 r 

pendix |A|) . The scaling T ~ Y 2 is unusual because it is v Jo v / 



in contradiction with the scaling relation (18), that was 
obtained with the use of Markovian scaling arguments. 
Having established the scaling law ( 65 ) in the Marko- 



V_e 2* — Yoe 2* 



7r 2 q 2 J °° dr ^-h-^e-^/W 



(68) 



vian approximation, we focus on the non-Markovian the- 
ory. Estimating the dependance of M q and Y n (r) with Yq 
for small Yo is not trivial: since there is an infinite num- 
ber of modes M qi the convergence of M q and YJ r (r) to 
as Yq — >• can be non-uniform. Indeed, the numerical 



A careful evaluation of these integrals using ( 66 ) leads to 



the corresponding form for the moments M q , valid under 
the hypothesis that < ft < 5/4: 



M q = g q 



(if 9 «i/y 2 ) 



(69) 
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where the coefficients g q are related to g{r) by: 
_ y/2 J^dr^- 3 / 2 (l-e-^)g(r) 



7r 2 q 2 



/ °° dr <Z,-i/2 e -i 2 * 2 



(70) 



For larger values of q, we define the new variable q = qY^. 
The behavior of M q then depends on the value of j3. Let 
us first assume that j3 < 1/4. In this case, we obtain 
M q ~ F 2+4/3 ft (5), with the function h(q) defined by: 



h(q) 



/2iJ~dTg(T)/MT)] 



3/2 



rig 2 | °° du u -i/4 e -g a ^ u -»(u)V(2V5) 



(71) 



Applying Eq. ( |62| ) leads to: 

/>OC 

Y w (t = uY o 4 ) = Y* dqh(q)(l 
Jo 



4/3 



(72) 



This expression is in contradiction with our initial as- 
sumption ( [66] ) : the case ^ < 1/4 is therefore not consis- 
tent with the theory. We now focus on the opposite case 
(3 > 1/4, for which we obtain that M q ~ F 2+4/3 h(q), 
with the function ft(g) given by: 



h(q) 



s/2 1 °° du u~ 3 / 4 (l - e~ q 77 ") e ^y(u) - e ^ 



TT 2 K f | °° du U -l/4 e -? 2 ^«-y(«)V(2^) 



(73) 



We note that the divergence g(r — » 0) ~ A/r^ is trans- 
ferred to the asymptotic form of g q (for large g) and h q 
(for small <?), as we have in these limits g q ~ g 2/3_1 and 



ftq ~ g 2/3 x . Using the relation (62), we obtain the fol- 
lowing links between g(r),y(u) and g q ,h(q): 

rOO 

y(u) = -y/2 / dq ft® (1 - e~^ u ) (74) 



s(r) - s(+oo) = V2^g q e-^ (75) 

Because g q ~ g 2/3_1 for large q, the series ^ ^ g is always 
divergent and leads to a divergent behavior of g(r) for 
small r. We identify this divergence as: 



g(r ->• 0) 



4^r(3/4) t/ 3 



(76) 



Initially, we had assumed that g{r) ~ A/r^, which 
is compatible with the asymptotic form ( |76] ) only if 
r(/3)r(l/4 - /3) = -4 v /7rr(3/4). It turns out that this 
equation has only one solution, that is j3 = 1/2. The 
value P = 1/2 is therefore the only value of j3 that is 
compatible with the theory. The non-Markovian theory 
in the limit Yq —> is completely defined by the cou- 
pled equations (70), (73), (74), (75) that are written in a 



consistent form that does not depend on Yq. Because 



f3 = 1/2 > 1/4, the equation for the reaction time (67) 
can be simplified, as only the part coming from the short 
time scales contributes: 



VN 3 / 2 ~ 



du 



[yW 



e 



(77) 



This expression is the most important result of this sec- 
tion, as it clearly shows that the reaction time scales with 



r3 

L 0_L_ 



contradiction with the Markovian approximation (165 



which predicts T ~ Yq. Note that the scaling T ~ Yq. is 
the scaling that is guessed by using simple (Markovian) 
arguments [see Eq. (18)]. The scaling (77) is supported 
by the numerical solution of the equations presented on 
Fig. [7[ This figure alone does not suffice to identify 
the limiting asymptotic behavior of T(Yq) because of the 
limited range of N where the numerical solution is avail- 
able (N < 320). For smaller values of N (N < 40), 
we had found that the non-Markovian theory is in close 
agreement with the results simulations, and it is there- 
fore likely that the non-Markovian asymptotic relation 
(77) is correct. Note however that is was derived in the 



framework of the stationary covariance approximation, 
and we do not know if the release of this approximation 
would change the scaling behavior of T. The calculation 
in this case is expected to be very cumbersome. 



E. The reactive shape of the polymer 

As stated above, the average shape of the polymer at 
the instant of the reaction is a key quantity that deter- 
mines the reaction kinetics. In this section, we give some 
information about what is the shape of the polymer at 
the instant of the reaction, especially when TV is large, 
and in the framework of the stationary covariance ap- 
proximation. The values of the moments M q are shown 
on Fig. [9] for a particular value of Yq = 2 and several 
values of N. It can be observed that when N becomes 
large, the moments M q reach an asymptotic curve which 
behaves as a power-law of q for large q. This power- 
law behavior breaks down when q/N becomes of order 
1, where finite size effects matter. The exponent of the 
power-law behavior M q that appears for large q can in 
fact be predicted by the theory. We show in appendix 
|E] that the only power-law that is consistent with the 
non-Markovian theory is: 



M q ~ -Moo/*/ 3 / 2 



(78) 



where is an unknown positive coefficient. This pre- 
diction is in agreement with the behavior of M q that is 
observed on Fig. [9] 

The knowledge of the average moments M q can be used 
to determine the average positions of the monomers at 
the instant of reaction (xi) n . Let us call Si the standard 
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s. ( = i/N ) 



FIG. 9. (color online) Average values of the modes M q (= 
rriq +1 /N) at the instant of reaction for Y = 2 in ID calculated 
in the stationary covariance approximation for various values 
of N. The oblique black line represents the scaling M q ~ 
^-3/2 rpj^ reac tive monomer is the first monomer. 



deviation of the position of the i th monomer at the in- 
stant of reaction. The theory predicts that there is a 
68% probability that the i th monomer is observed be- 
tween (xi)^ — Si and (xi)^ + Si. These two curves are 
represented on Fig. [lOj together with the average reac- 
tive shape of the polymer and two examples of polymer 
reactive conformations. As can been observed on this 
figure, the non-Mar kovian theory predicts a significative 
shift with respect to the reactive point of the positions of 
all the monomers at the instant of reaction the reactive 
non-equilibrium conformations of the polymer is very dif- 
ferent from an equilibrium conformation (for which the 
average positions vanish: (xi) n = 0). It can also be seen 
on this graph that the curve (x(s)) 7T (where Si is the 
position of the i th monomer in the chain) shows sharp 
variations for small values of s. The origin of this anoma- 
lous behavior of (x(s)) 7T is due to the power law behavior 
(78) of the co effic ients M q . Indeed, taking the continuous 
limit of Eq. (|5|6|), we obtain the relation: 



N 



-y/2^M q [l -cos(sTrg)] 



(79) 



q=l 



Inserting the asymptotic behavior M q ~ 1/q 3 / 2 and re- 
placing the sum by an integral yields, for small s: 



N 



V2M C 



f 

Jo 



dy 1 — cos(?/7r) 
s (y/s) 3 / 2 



(80) 



The asymptotic behavior ( 80 ) indicates that the slope of 
(x(s)) 7r is infinite at the point 8 = 0: the reactive position 
of the first monomers of the chain is therefore strongly 
shifted with respect to the position of the reactive region. 

Finally, we describe the dependance of the reactive 
shape of the polymer on the initial distance between the 



FIG. 10. (color online) Positions of the monomers at the in- 
stant of reaction in ID. We show the average reactive position 
of the monomers (continuous red line) (x(s)) 7T as a function 
of the coordinate s = i/N of the monomers in the chain. 
The dashed red lines are (x(s)) 7r ± S(s), with S(s) the stan- 
dard deviation of x(s) predicted by the non-Markovian the- 
ory. We also show two examples of reactive conformations 
generated from the theoretical distribution of reactive con- 
formations (green and blue curves). Parameters: N = 320, 
Yb = 2 and p = 1. 



reactants. On Fig. 11, we represented the average po- 
sitions of the monomers for several values of Yq. As Yq 
is decreased, we observe the apparition of two regions in 
the curve (x(s)) 7r . There is a small region around 8 = 0, 
whose size decreases with Yq, in which (x(s)) 7r does not 
depend much on Yq. There is another region, for larger 
values of s, where (x(s)) 7T varies slowly with s but de- 
pends strongly on Yq. The presence of these two distinct 
regions is the sign that distinct length and time scales 
appear in the problem when Yq — » 0, and is in agreement 
with the structure of boundary layer ( 66 ) that was pos- 



tulated in section llV Dl A similar structure is observed 
for the coefficients M q : the large q part of the spectrum 

tis independent on Yq, but disappears as Yq — » (Fig. 
inset). The fact that the small length scales part of 
s))^ and the large q part of M q is independent on Yq 
can be predicted by the theoretical analysis. We have al- 
ready seen that for small Yq and large g, the correct scal- 
ing law for M q is M q c± Y 3 h(qY 2 ). To be consistent with 
the scaling law (78), we must have h(qY 2 ) ~ l/(qY 2 ) 3 / 2 , 



from which we deduce that M q ~ l/q 3 / 2 : M q is asymp- 
totically independent on Yq for large q. 



F. Concluding remarks on the ID problem 

At this stage, we have exposed a description as com- 
plete as possible of a non-Markovian theory that enables 
the determination of the first passage time of a monomer 
of a Rouse polymer chain to a given target in a one 
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FIG. 11. (color online) Average reactive position of the 
monomers in ID for N — 320 and several values of Yq, when 
the reactive monomer is the first one. In the inset, we show 
the coefficients — M q for the same parameters (the color code 
for the curves is the same). 



dimensional space. The key approximation of the non- 
Markovian theory is that the distribution of the polymer 
conformations at the instant of reaction is a multivariate 
Gaussian. The non-Markovian theory and its simplified 
version (that uses the stationary covariance approxima- 
tion) are in good agreement with numerical simulations, 
to the difference of the Markovian approximation, which 
assumes that the polymer is at equilibrium when the re- 
action occurs. One of our most important results is that 
the Markovian approximation predicts a different asymp- 
totic relation for the reaction time T with the initial dis- 
tance Xo. As the non-Markovian theory is supported by 
simulations for the values of N that we tried, we deduce 
that the Markovian approximation predicts reaction time 
that are very largely overestimated for long chains. We 
have also described the asymptotic behavior of the aver- 
age reactive polymer conformations, and we have showed 
that their spectrum are characterized by a slowly decreas- 
ing power-law tail. 



V. NON-MARKOVIAN REACTION KINETICS 
IN 3D 

A. Generalization of the theory to 3 dimensions 

Up to now, we have only considered the case of a one- 
dimensional space. However, the theory can be extended 
to the case of a d— dimensional space. We now describe 
the non-Markovian theory for a 3-dimensional space, but 
we will give less details than for the theory in ID. Note 
that the theory for a two-dimensional space, for example, 
could be easily obtained by following the successive steps 
of our approach. There are some differences between the 
ID and 3D situations that must be taken into account to 
properly write a non-Markovian theory in 3D. The first 



difference with the ID case is that the size of the spheri- 
cal reactive zone has now a finite radius a, and therefore 
we have to consider the "entrance direction" that is de- 
fined by the direction of the vector R at the instant of 
reaction. This direction defines the azimuthal angle 6 
and the polar angle cp at the instant of the reaction, and 
we use the notation Q as a shortcut to represent y?, 6. A 
second difference with the ID case lies in the choice of 
initial conditions. In the case of a large confining vol- 
ume, we anticipate that the reaction time depends only 
on the distance Rq between the reactants. Therefore, it is 
not restrictive to assume that the initial configuration is 
isotropic. Specifically, we assume that initially the poly- 
mer is at stationary state, with the restriction that the 
distance between the reactive monomer and the center of 
the target is Ro. In this case, the initial distribution is 
a superposition of Gaussian distribution, averaged over 
angles: 



|a)) = J dfiP 8tat (|a>|i?ou r (n)) 



(81) 



where dVt = sin QdQdtp/ (An), and u r (Q) denotes the ra- 
dial unit vector pointing outwards the reactive sphere. 
The third difference with the ID case resides in the way 
of writing the renewal equation. Let us consider a poly- 
mer that is observed at t in a configuration |a) with the 
reactive monomer at position Rj ((b|a) = R/). The pa- 
rameter Hf can be chosen arbitrarily inside the reactive 
region. Observing the conformation |a) at time t nec- 
essarily implies that the polymer has reached the target 
for the first time at some time t' <t, with some entrance 
angle Vt = (9,<p) and some configuration |a') (that is such 
that (b | a') = au r ). Therefore, if we define fo(\a!),t') as 
the probability density that the reactive region is reached 
for the first time at t' with a configuration |a') given that 
the entrance angle is Q! , we can write the following re- 
newal equation: 

P(|a),t|{ini},0) = 

f dt' J d\s!) J dto' fo(\s!) J)P(\*),t-lf\ |a'),0). 

(82) 

This equation takes into account the fact that the reac- 
tion can occur with equal probability at any place on the 
reactive sphere. We introduce the probability 7r^(|a)) 
of reacting with a configuration |a) given that the re- 
active monomer position R has the angular coordinates 
Q = (0, if) when the reaction takes place. As in the ID 
case, taking the (temporal) Laplace transform of (82), 



expanding for small values of the Laplace variable and 
taking into account the superposition relation (81) leads 
to: 

rP st at(|a>|R/) (Rt) = 

dt J dQ [P(|a),t|7r n ,0)-P(|a),t|{i? u r ,8tat},0)] 

(83) 
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As in the ID case, we make a large volume approxima- 
tion: all the terms that are appearing in Eq. (83) are 



approximated by their value in infinite space, except for 
the term P stat (Tlf) = 1/V. We also make a Gaussian 
approximation of the splitting probability distribution 
ttq. Writing a complete theory requires to determine 
a 3 N x 3N covariance matrix and a 3 N mean vector. 
Here, for simplicity, we restrict ourselves to the simple 
case where the covariance matrix of each spatial coordi- 
nate of ttq is given by its stationary value: 



coy (a iia aj^) 



r stat,* 



(84) 



where a, j3 stand for spatial coordinates. This approxi- 
mation is the equivalent to the "stationary covariance ap- 
proximation" that was developed in the ID case. In the 
approximation ( 84 ) , we have assumed that the covariance 



matrix is isotropic. Now, for symmetry reasons, only the 
radial components of ttq can have a non- vanishing mean 
vector, so that we can define the average radial modes at 
the reaction mf with the relation: 



(85) 



The self-consistent equations that define the values of 
are obtained by multiplying Eq. (83) by a iz and by 



integrating over all the modes. The detailed calculation 
is presented in the appendix [FJ and we arrive at the fol- 
lowing equation, valid for i > 2: 



dti 



R^mje 



R%\ bjjl-e-W 
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(86) 



Note that the equation that defines the moments mj de- 
pends in general on the choice of the parameter Rf, which 
can be arbitrarily set between and a. The expression 
(86) corresponds to the particular choice Rf — >• 0. In 



Eq. (86), we have used the notation R n (t) to represent 
the average "radial" position of the reactive monomer at 
a time t after the reaction. R n (t) is in fact the average 
position of the reactive monomer in the direction defined 
by the entrance angle to the reactive region, at a time t 
after it has reached the reactive zone for the first time. 
It is given by: 



N 



R^(t) = {b\^)=a-Y J ^<^ 



-Xit\ 



(87) 



i=2 



The reaction time reads: 



T 
V 



J dt J dTt\ exp 



(Rf - R^WrY 



exp 



(Rf - R u r ) 2 

2i/j 



Because of isotropy, we can assume without loss of gen- 
erality that R/ is located on the z-axis: R/ = RfU z . 
Noting that (R f - R n u r ) 2 = R\ + (R f ) 2 - 2R 7T R f cos0, 



we can integrate (88) over the angles: 
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This expression is simplified by taking Rj = 0: 
dt 
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(89) 



. (90) 



The set of TV — 1 self-consistent equations ( 86 ) together 
with the expressions of the reaction time (89), (90) com- 



pletely define the non-Markovian theory in 3D under the 
stationary covariance hypothesis. As in the ID case, we 
can also define a Markovian approximation, where the 
splitting distribution is approximated by an equilibrium 
distribution: 



TIL; 



Sna/bi ]R n (t) 



(Markovian Approx.) (91) 



|90 ) gives the 



Reporting this approximation into Eqs. 
Markovian expressions for the reaction time. Note that in 
fact, both Markovian and non-Markovian theories do not 
predict a single value of the reaction time, as the result 
depends on the parameter Rf, which can be in princi- 
ple chosen arbitrarily with the restriction < Rf < a. 
The fact that the final result depends or not on Rf is a 
test of consistence of the theory. We will see below that 
the Markovian approximation fails to pass this test, as 
it gives two different values corresponding to Rf = or 
Rf = a, and these two values are not upper and lower 
bounds of the correct result. For the non-Markovian the- 
ory, the numerical integration of the equations shows 
that the value of the mean first passage time obtained 
for different Rf are almost undistinguishable. The non- 
Markovian theory is therefore consistent. 



B. Comparison with numerical simulations 

In order to characterize the validity of the non- 
Markovian theory, we compared its predictions with the 
results of stochastic simulations. On Fig. [l2j we repre- 
sented the average reaction time for a moderate value 
of TV. As we can see on this figure, there is a very 
good agreement between the non-Markovian theory and 
the simulations. The Markovian approximation is qual- 
itatively correct, but does not match quantitatively the 
data. In addition, it gives two distinct results correspond- 
ing to the two possible choices of Rf. As in the ID case, 
renormalizing the values of ai at the instant of reactions 
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by their means and variance and superposing all the his- 
tograms gives a curve that is very similar to a Gaus- 
sian function 



(Fig. 



15), suggesting that the Gaussian 
approximation is a very accurate one. However, the non- 
Markovian theory does not predict the correct values of 
the position of the polymer center-of-mass at the instant 



of reaction rril/yN (Fig. 13). This discrepancy could 



possibly come from the stationary covariance approxi- 
mation (84), which assumes in particular the isotropy of 
the covariance matrix. In fact, the coefficient cr^/N is 
very well approximated by its stationary value in the ra- 
dial direction, but is underestimated by about 25% in the 
perpendicular directions (see Fig. [l4| . We conclude that 
the non-Markovian theory (with the stationary covari- 
ance approximation) provides an accurate description of 
the reaction time in 3D, although there is a disagreement 
between the theoretical and measured values of the mo- 
ments mf, afj. 



0.35 

0.3 

0.25 

> 0.2 
I- 

0.15 
0.1 
0.05 





o Simulations 
— non-Markovian 
Markovian(R f =a) 

Markovian(R f =0) 



10 



20 



30 



40 



FIG. 12. (color online) Comparison between the values of the 
reaction time predicted by the theory and measured in the 
simulations in 3D for N = 20. Red circles: simulation results 
(the error bars are 95% confidence intervals ; each symbol 
is the result of an average over 1470 simulation runs). The 
capture radius is a = 2, and is located in the center of a 
spherical confinement volume of radius R = 40.25. The time 
step is At = 0.0005. The two green upper curves (dashed and 
continuous) represent the Markovian estimates of the reaction 
time and correspond to the choices Rf = and Rf — a. 
The red continuous curve represents the results of the non- 
Markovian theory (for which the two curves for Rf — and 
Rf — a are almost superposed). 



C. The polymer reactive conformations 
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FIG. 13. Position of the polymer center-of-mass at the instant 
of the reaction in 3D for N = 20 in the parallel direction (the 
direction defined by the entrance angle to the reactive region) . 
Red circles: results of simulations. Red line: non-Markovian 
theory. Dashed line: Markovian approximation (in which the 
center-of-mass is located at the surface of the target, at a 
radial position a = 2). All parameters are the same as in Fig. 
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FIG. 14. (color online) Variance of the polymer center-of- 
mass position at the instant of reaction for N = 20 in 3D, in 
the radial directions (in which the reaction takes place, lower 
curves with circles) and the 2 perpendicular directions (upper 
curves with squares and diamonds) . The dashed green line is 
the prediction of both Markovian and non-Markovian theories 
in the stationary covariance approximation. The difference 
between these two directions is a signature of the presence 
of a weak anisotropy in the covariance matrix, which is not 
accounted for in the theory. All parameters are the same as 
in Fig. [12] 



1. Limit of small target size 

We now study the solution of the equations of the non- 
Markovian theory in various limiting cases. For simplic- 
ity, we restrict ourselves to the case where the initial 



distance between the reactants is large (Rq — ^ oo), in the 
regime where T does not depend on Ro anymore. We 
first focus on the case a — » for a fixed value of N. Let 
us assume that the moments mf vanish when a — » 0, and 
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a. (centered and rescaled) 



FIG. 15. (color online) Superposition of all the histograms 
of ai (renormalized by their means and variance) in all the 
directions and for all values of Ro, superposed with the stan- 
dard Gaussian distribution. The data are the same as in Fig. 
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that they are proportional to a. Then the function R^(t) 
is also proportional to a and the simplification R n (t) ~ a 
is correct at short times t —> 0. In the limit a — >• 0, all the 
integrals appearing in Eq. (86) are dominated by their 



short time part: they can be estimated by approximating 
the integrands by their short time limit. For example, us- 
ing the simplifications R^{t) — a, ip — 2t and e~ Xit ~ 1, 



we get for the first term of Eq. ( 86 ) : 
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(92) 



Using the same simplifications (and 1 
evaluate the second term of Eq. (I86|) . 
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(93) 



By Eq. (86), the expressions Eqs. (92) and ([93]) must 
compensate each other, and we obtain ~ — abi for 
i > 2. From the condition (blm^) = a, we also get the 
result ml = a\/N(2 — l/N). The fact that the moments 
are proportional to a validates our analysis. From this 
expression, we deduce that the average radial position 
are (zi) n = a(l + #i,i): the monomers are therefore not 
located at the surface of the reactive zone on average. 
Because the moments are proportional to a, they do 
not enter in the simplification of the expression of the 
reaction time at lowest order in a, which reads: 
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This expression is obtained by approximating the inte- 
grand of Eq. (90) by its short time limit. This result is 



valid for both Markovian and Non-Markovian theories in 
the limit a — >• 0, which is an indication that both theories 
predict the good result in this limit. It is also consistent 
with the scaling relation (16). Equation (|94|) shows that, 



for a very small size of the reactive region, the reaction is 
limited by the time that a single monomer, disconnected 
from the rest of the chain, finds the reactive zone. 



2. Limit of large N 

We now consider the limit of a large number of 
monomers: we assume that N — oo when the param- 
eter a = a/y/~N remains constant. As in the ID case, 
we have X q ~ (q — l) 2 tt 2 / TV 2 , b q ~ \/2/iV, and we intro- 
duce a rescaled time r = t/N 2 . We assume the scaling 



M q — rriq +1 /N, which is the scaling for which Eq. (86) 
does not depend on N any more, as it becomes: 
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(95) 



where ^ is given by ( 63 ) and the rescaled function (r) 
reads: 
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The evaluation of the mean first passage time is : 
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= F(a) 



(96) 



(97) 



Here, F is a dimensionless function that depends only 
on a = a/y/N (because Y n itself depends implicitly on 
a). In the Markovian approximation, where — a, the 
function F has a simple expression: 
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This function can be developed for small values of the 
reactive zone: 
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0.112 (a ->• 0) 



(99) 



Hence, in this regime, the reaction time does not depend 
on the size of the target: it is a a signature of the compact 
search of the monomer at short time scales. Mathemati- 
cally, it comes from the subdiffusive behavior of the mo- 
tion at short time scales that implies that ^ ~ t 1 / 4 , and 
therefore makes the integral (99) a convergent one. In- 



terestingly, the numerical coefficient (99) is the result of 
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an integration that runs over short and large time scales 
r, and therefore depends on the properties of the motion 
at all time scales. This remark is fully consistent with 
the analysis above Eq. (14), where we had found that 



the reaction time is the result of two substeps (one dif- 
fusive at large time scales, one subdiffusive at short time 
scales) that last approximately the same time. 

In 3D, the Markovian and the non-Mar kovian theories 



predict the same scaling law (97), but the dimensionless 
function F is different. In the non-Mar kovian theory, F 
has to be determined numerically, and it is represented 
on Fig. 16 The asymptotic behavior of M q for large g, 



however, can be analytically determined. Let us consider 
Eq. (95) as G q ({Mi}) = 0. We consider the development 
of G q in powers of q. The first term of this development 
is of order 1/g 2 , its coefficient must vanish, leading to: 
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This is a global relation that involves all the moments 
M q . Inserting this equality into Eq. (95) leads to the 



estimate of M q as a ratio of two integrals: 

q ttV Jo°° dTY^-^e-* 2 * 

with the function H defined by: 
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In the general case, the expression (101) is not sum 



cient to determine the M q , because Y n does depend on 
M q . However, due to the presence of the term e~ q2 ^ T ^ 
for la rge q. it is clear that the integrands appearing in 
Eq. (101) can be evaluated at their short time limit, 
where Y n ~ a. The integrals appearing in ( 101 ) can 
be calculated with the saddle point method. Solving for 
H'(t*) = 0, we get the position of the saddle point at 
r* ~ [a 2 /( 4<7 2 ^7 r 2 )1 2/3 in the limit q -» oo. Then, the 
expression (|101|) can be evaluated as: 
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(103) 



This result is in good agreement with the computed val- 
ues of M g , even for reasonable values of q and TV, as can 
be seen in Figure [l7|), where the values of M q differ from 
the asymptotics ( |103| ) by a factor smaller than 1.4 for 
2 < q < 400. The fact that the coefficients M q decrease 
as a slow power-law of q implies that the function (z(s)) 7r 
does not admit a derivative around s = 0. Mo re prec isely, 
using the same method as in ID [see Eqs. f79p0l )], we 
get: 
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(104) 



This formula means that the monomers that are close 
from the reactive monomer in the chain have a position 
at the instant of reaction that is significantly shifted with 
respect to the po sitio n of the reactive site. When a = 0, 
the scaling law (103) is not valid any more. Prelimi- 



nary analysis suggests that in this case the asymptotic 
behavior of M q is still characterized by a power-law, and 
becomes M q ~ (In q) 1 / 2 /q 3 / 2 . 
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FIG. 16. (color online) Reaction time in 3D for various values 
of N as a function of the rescaled capture radius a = a/y/N in 
the Markovian approximation (upper dashed line) and non- 
Markovian theory (lower curves with symbols) for various val- 
ues of N. The thick dark lines are the Markovian (dashed) 
and non-Markovian (continuous line) estimates of the scaling 
function F(a) that is reached for large N. The divergence of 
the rea ction time for small a is due to the asymptotic behav- 
ior (94). All Markovian estimates use the value Rf = 0. The 



reactive monomer is the first monomer. 



D. Effect of the monomer position in the chain 

Up to now, we have considered only the case where the 
reactive monomer is the first monomer (p = 1). However, 
the equations of the non-Markovian theory are written 
for any value of the position of the reactive monomer 
p (which enters in the definition of the coefficients bi). 
We now complete the study by briefly studying the ef- 
fect of the position of the reactive monomer in the chain. 
The reaction time as a function of p is represented on 
Fig. 18 in the case of a large initial distance between 
the react ants. As can be observed, varying the position 
of the monomer does not have a dramatic effect on the 
reaction time, but it is clear that the reaction time is 
reduced when the reactive monomer is located close to 
the polymer extremities. This observation can be un- 
derstood by considering that the motion of an exterior 
monomer is less hindered by the polymer chain in the 
subdiffusive regime, as they are surrounded by only one 
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FIG. 17. (color online 
reactive region a = a/y/N = 0.1. 
right correspond the the values N 
thick dashed black line represents the expression (103) and 
has slope —4/3. The reactive monomer is the first monomer. 



Coefficients M q in 3D for a size of 
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FIG. 19. (color online) Average positions of the monomers 
at the reaction in 3D, for 3 different values of the index of 
the reactive monomer (p = 1, p = 41, p = 101). Parameters: 
N 200, a 2.8284 and R oo. 



VI. CONCLUSION 



polymer chain (instead of two chains that are surround- 
ing the interior monomers). This faster motion at small 
time scales leads to a smaller reaction time. The differ- 
ence between the results of the Markovian approxima- 
tion and the non- Markovian theory is maintained when 
the reactive monomer is moved along the chain. We also 
represented the polymer reactive shapes for different val- 
ues of p on Fig. [l9] one can observe that the shape (xi)^ 
has a singular behavior around p = z, a fact which is re- 
lated to the slowly behavior of the coefficients ra£ as a 
power- law of q. 
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FIG. 18. (color online) Reaction time in 3D for several values 
of the index of the reactive monomer in the chain p. Param- 
eters: N = 200, a = 2.8284 and R -> 00. 



In this paper, we have presented a theory that de- 
scribes the kinetics of intermolecular polymer reactions 
in the diffusion controlled regime. The theory takes ex- 
plicitly into account the non- Markovian nature of the 
monomer motion by determining the distribution of the 
polymer conformations at the very instant of the reac- 
tion. The key hypothesis of the theory is that this dis- 
tribution is a multivariate Gaussian, which enables the 
derivation of a set of self-consistent equations that de- 
fine the parameters of the distribution of reactive con- 
formations. Another hypothesis of the theory is the 
large volume approximation, and our study generalizes 
approaches that use this approximation in the case of 
Markovian processes [36 j. Comparison with the results 
of numerical stochastic simulations shows that the non- 
Markovian theory predicts very accurately the reaction 
time, both in one dimensional and three dimensional 
spaces, and for all the values of parameters of the prob- 
lem (number of monomers, size of the reactive region, 
initial distance between the reactants and position of the 
reactive monomer in the chain). The non-Mar kovian the- 
ory gives much more precise results than the Markovian 
approximation, in which the distribution of reactive con- 
formations is replaced by the polymer equilibrium dis- 
tribution. This Markovian approximation is equivalent 
to the Wilemski-Fixman approximation in the context of 
intramolecular reactions [8[ |9j [251 ED] , and is also similar 
to the approximation of quasi-independent intervals [28] 
in the context of general Gaussian processes. The distri- 
butions of reactive conformations predicted by the non- 
Markovian theory are in general very close from the ones 
measured in simulations, and it is in fact quite surprising 
that the marginal laws for the Rouse modes at the instant 
of reaction are very close from a normal distribution. 
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We have also described a simplified non-Markovian the- 
ory, the "stationary covariance approximation", which 
catches the main non-Markovian effects and is in close 
agreement with simulations. 

In addition, we have derived various asymptotic be- 
haviors of the reaction time in the non-Markovian and 
Markovian theories. One of the most interesting result 
of our study is that it reveals a strongly non-Markovian 
regime in ID, where the Markovian theory predicts a 
wrong asymptotic relation of the reaction time as a 
function of the initial distance. In this regime, the re- 
active conformations are so different from the equilib- 
rium conformations that the Markovian approximation 
leads to expressions of the reaction time that can be 
overestimated by several orders of magnitude for long 
chains. When the reaction occurs in 3D with long chains, 
the non-Markovian effects are only quantitative and the 
Markovian theory overestimates the reaction time by 
roughly 30% - 100%. 

We have also shown that one can derive scaling rela- 
tions for the reaction time in a systematic way by con- 
sidering the decomposition of the reaction into different 
substeps that occur at several length and time scales, 
where the properties of monomer dynamics are different. 
Despite the fact that these guesses are based on Marko- 
vian arguments, the scaling relations obtained in this way 
are always in agreement with the non-Markovian theory. 
These reasonings can help to identify which substeps of 
the reaction involve diffusive or subdiffusive regimes, and 
which are the dominant substeps. Interestingly, we no- 
tice that, when the dominant substeps are diffusive, the 
Markovian and non-Markovian theories predict the same 
value for the reaction time at lowest order, a fact which 
must be closely related to the fact, among non- stationary 
Gaussian processes, only Brownian motion is Markovia. 
Each time the scaling arguments predict that at least 
one of the substeps is subdiffusive, the difference between 
Markovian and non-Markovian theories is at least quan- 
titative: it is the case in 3D, when the capture radius is 
small compared to the polymer size, but not that small 
so that the polymer can be considered as a continuous 
chain. As stated above, in ID the situation is more ex- 
treme, as the Markovian and non-Markovian expressions 
for the reaction time can differ by several orders of mag- 
nitudes. 

In this study, we have also described what is the typical 
shape of the polymer at the instant of reaction. The reac- 
tive polymer conformations are much more elongated on 
average than equilibrium conformations. Because of this 
elongation, the reactive monomer explores more space 
around the polymer center-of-mass than in an equilib- 
rium configuration, leading to a faster reaction kinetics. 
This picture also holds when the reactive monomer is in 
the interior of the chain. The kinetics is slowed down 
when the monomers are in the chain interior, due to the 
fact that the motion of an interior monomer is hindered 
by two polymer chains instead of only one for exterior 
monomers. Our analysis reveals that the reactive con- 



formations are characterized on average by a slowly de- 
creasing tail in the spectrum: the average values of the 
Rouse modes (or, equivalently, the Fourier coefficients of 
the average polymer shape) decrease as a power-law of 
the wave number with an exponent that can be calculated 
analytically. As a consequence, the average position of 
the monomers at the reaction shows a singularity around 
the reactive monomer, meaning that the monomers that 
neighbor the reactive monomer are significantly shifted 
from it at the instant of reaction. 

In this paper, we have investigated only the case of 
intermolecular reactions involving Rouse chains. Due to 
the multiciplicity of the time scales involved and the non- 
Markovian nature of the problem, determining the pre- 
cise reaction kinetics in this case is not a trivial task. We 
hope that our non-Markovian approach can be applied 
in the future to understand the non-Markovian effects on 
the kinetics of reactions that involve more complex poly- 
mers such as branched polymers [43 or polymers with 
excluded volume [44 or hydrodynamic interactions [45] . 
It would also be interesting to generalize the theory to 
investigate first passage time properties in the more gen- 
eral context of non-Markovian Gaussian processes. 
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Appendix A: The function ip(t) for large N. 

In this appendix, we describe one way of obtain- 
ing the asymptotic behavior [Eq.([l2])] of the function 
ip(t) for large values of TV and intermediate time/length 
scales. First, we have to distinguish between a a reactive 
monomer located at the interior or the exterior of the 
chain. For this sake, we introduce the parameter s that 
describes the position of the reactive monomer when one 
takes the limit of infinite N: 



lim 



P-1/2 
N 



(Al) 



If s = or s = 1, we say that the reactive monomer is 
located at the exterior of the chain (it is surrounded by 
only one polymer chain). If < s < 1, the monomer 
is at the interior. For any value of s, ip(t) reaches an 
asymptotic form for large N: 



(A2) 



where ^ is a function that depends on the rescaled time 
r = t/N 2 that can be easily identified with Eqs. (p5|,(10): 
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(A3) 
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The function \I/ obviously vanishes at r = 0, but does 
not admit a finite derivative in r = (be caus e the inser- 
tion of (1 - e- q2 ^ T ) ~ q 2 7T 2 r into Eq. |A3| leads to a 
diverging series). In the limit r — ^ 0, if < s < 1, the 
term [cos(sqTr)] 2 varies very fast compared to the other 
terms of the series, and it can therefore be replaced by its 
average value 1/2. If 8 = of 5 = 1, the term [cos(sg7r)] 2 
is trivially replaced by 1. When r — >• 0, the variable 
V = qVt can be considere d a s a continuous variable. 
Then, replacing the series (A3) by an integral, we ob- 
tain: 




if < s < 1 
if s = 0, 1 



(A4) 



A simple evaluation of this integral leads to the behavior 
\I/ ~ /it for interior monomer and \I/ c± ^\Jr Jtt for an 
exterior monomer. Rescaling \I/ and r by the appropriate 
powers of N leads to the expression ( 12 ) in the main text. 



Appendix B: Projection formulas 

Here, we briefl y describe how to derive the projection 
formulas ( 42|43 ), which are an adaptation of a result on 
conditional gaussian distributions that can be found for 
example in the chapter 3 of Ref . [41] . Consider a set of 
N gaussian random variables ai, a/v, with mean vector 
rrii and covariance matrix 0^-. We consider the average 
of ai given that another variable ctj takes the value a®. 
According to Ref. [41], 



6i,...,6jv- If b n is the only non- vanishing coefficient, it 
is trivial that these formulas are also true for the A^ th 
variable. If it is not the case, doing the same reasoning 
with another variable that has a non- vanishing coefficient 
leads to the conclusion that the formula is also true for 
the N th variable. 



Appendix C: Numerical solutions of the equations of 
the non-Markovian theory 

Obtaining theoretical estimates of the non-Markovian 
theories requires to be able to solve a system the 
system of equations ( 32|33 ). Let us write Eq. 
Q as Gi({m]},{aJ k ^i = and Eq. Q as 
j'Iij({rrik} : {cr kl }) = 0. In order to obtain the solutions of 
these equations, we introduce a fictive "time" s and we 
numerically solved the dynamical system: 



d 
ds 
d 



ds^ 



mr = -G,({mJ},KJ) 
-^({<},K/}) 



(CI) 
(C2) 



The solution of this dynamical system converges to the 
solution of the non-Markovian theory as s — )• oo, if one 
takes an initial solution that is such that (b\m n ) = and 
a^lb) = |0). In practice, we take the values of mf at,X ° 
and a ^ at > x o a s initial conditions. Another difficulty is the 
presence of indefinite integrals in the equations ( 32|33 ). 
These integrals are numerically evaluated in Matlab with 
a vectorized adaptive Gauss-Kronrod quadrature algo- 
rithm [46] . 



E(a i |a fe = A) =rrii — (m k - A) 



%k 



(Bl) 



Similarly, the formula for the covariance of a^o^ given 
that the variable a k takes the value A reads: 



cov(ai,a,j\ak = A) = k 



Okk 



(B2) 



Now, consider the variable X = (b\a) such that at least 
one of the bi is different from 0. Assume that b n ^ 0. 
Then, consider the other distribution P(ai, ajv-i, X). 
This distribution is also gaussian, the average of X being 
(b\m), while the covariances read: cov(<2i,X) = (b\0)ei 
and var(X) = (b\0\b). Then, we apply the two formulas 
( |B1[ ),(B2) to the variables (ai, ajv-i, X) (we replace 
m k by~Jb\m), kk by var(X) and ik by cov(a i ,X)): 



E(a z \X = X ) = m % - ((b\m) - X ) 
cov(a i ,a J |a fc = A) 



(ei\6\b) 



(B3) 



(b\e\b) 

{e l \e\b){e J \e\b) (B4) 



(b\e\b) 

These relations are exactly the projection formulas 
( [42143] ) of the main text: they are true at least for 



Appendix D: Simulation algorithm 

In this appendix, we describe the method of simula- 
tions that we used in ID. In a simulation run, the initial 
value of each mode (i > 2) is taken from a normal 
distribution of variance 1/A^. The initial positions of the 
monomers are simply obtained by applying Eq. (|6| and 
by translating the whole polymer to that the initial po- 
sition of the first monomer is Xq. Then, at each time 
steps, when the polymer is far from the absorbing and 
the reflecting wall, the positions evolve according to: 



Xi(t + At) =Xi(t) 



N 



(Dl) 



in which u\n, un are N random number taken from 
a centered gaussian distribution with variance 1. When 
the reactive monomer is close t o th e absorbing wall or 
to the reflecting wall, equation (Dl) is not very precise, 



as stated by Peters et al [42] . In the case of the proxim- 
ity with the reflecting wall, it misses the fact that there 
is a shift towards outside the reflecting wall because it 
does not take into account the fact that the particle has 
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a decreased probability to approach the wall (and zero 
probability to cro ss it ). Following Peters et al [42], we 
then modify Eq. (Dl) into: 



The slowest term is of order q 2 and must vanish, which 
implies the global condition: 



x 1 {t+At)= x^-fl 



refl / L ~ X l 



At 



f 

JO 



dr 

^3/2 



exp 



Y 2 ' 

-21 1 ,; 



exp 



I Y 
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"I /: 



refl / L ~ X l 



= 

(E3) 



At k(x\ X2)At (D2) Inserting this relation into Eq. (IgTI) , and rearranging the 



In this equation, u\ is a random number that takes the 
values ±1 with equal probability, and the positive func- 
tions /{ efl and f™^ are the functions fi and f 2 of the 
equation (18) in the reference [42 j. The supplementary 
terms take into account a shift in the direction opposite 
to the wall. Note however that they have been calculated 
by explicitly solving the Fokker-Planck equation near a 
reflecting wall in the case of a single particle, which is not 
the case here because of the presence of many monomers. 



However, we still expect that (D2) is a good approxima- 
tion of the dynamics in the limit At — >• 0. 

When the first monomer is close from the absorbing 
wall, one first calculates P a b s = 1 — eri(x\ / {2\f~At)) the 
probability of being absorbed between t and t + At. One 
then generates a random number between and 1 to 
decide whether or not the target is reached during the 
time step, in which case the simulation stops. If the 
absorbing wall is not reached, then x\ evolves according 
to: 



Xl (t+At) = Xl (t) + fl 



abs / X l 



At 



Ul /f bs ( -p= ) At - kfa - x 2 )At (D3) 



where the random number u\ takes again the values ±1 
with equal probability. The positive functions /f bs and 
yabs are £ ne f unc tions f\ and f 2 of the equation (16) 
in the reference [42]. At the end of the simulation, the 
positions of the monomers are recorded, thereby giving 
an access to the splitting probability. 



Appendix E: Asymptotic behavior of M q in ID 



In this appendix, we prove that the only power-law 
behavior of M q that is compatible with the theory is 



M„ 



-3/2 



Let us postulate the form M q ~ —M^/q 1 
oo. The fact that M q is a summable series im- 



We also assume that 7 < 3. Replacing 



for q 

poses 7 > 1 

the sum by an integral in the expression ( 62 ) , we get the 
short time behavior of Y n (r) : 



Y n (r) — t— >o A-T a ; a=( 7 -l)/2 



(El) 



A=-V2 dy(l 
Jo 



e~y 2 - 2 )/y-y (E2) 



Let us consider Eq. (61). All the terms of the develop- 



ment of its right hand side in powers of q must vanish. 



remaining terms leads to the formula: 

y/2 / °°dr m-le- q2 ^ T (V^e-^f -Y Q e~^ 



M q = 



7r 2 q 2 J^dr ^-f e -^ 2 -e-^ 2 /(2*) 



(E4) 



First, we use the saddle point method to evaluated the 
contribution of the integral that depends on Yq. We 
write: 



poo Y 2 rc 

Jo Jo 



dr ^"2 e - 
where the function G is: 
G = 7r 2 r + 



dre 



- q 2 G 



Y 2 
x o 



2q 2 K^ 



31nr/4 



(E5) 



(E6) 



Solving for G'(r*) = yields the value r* = 
[Yo 2 /(4(7 2 ^7r 2 )] 4 / 3 - 1/g 2 / 3 (for large q). Therefore, the 
saddle point method implies that: 



f 

Jo 



dr \&-i e - q27T2r -™ - e-^G^n „ e~ q 



4/3 



(E7) 



To evaluate the other integrals of (E4), one must dis- 
tinguish between the cases a < 1/4 and a > 1/4. Let 
us first consider the case a < 1/4. Then, Y 2 j^l ~ r 2a - 2 
diverges for small r. We pose H(r) defined by: 



H(t) 



7T 2 T ■ 



Y 2 

7T 

2q 2 ^ 



: 7r 2 r ■ 



-T 



2a-l/2 



(E8) 



Solving for H'{r*) = leads to the following position of 
the saddle point: r* ~ 1/ '^ 1 /( 3_4a ). Hence, the saddle 
point method indicates that: 



f 

Jo 



T 3 
\I/2 



e~ q 71 r F 7r e"2f e 



2-4a 
7 3-4a 



(E9) 



Comparing with ( E7 ) , it is clear that this term dominates 
the term that depends on Yq in (E4). Omitting the term 
that depends on Yq in (E4) and using the saddle point 
method leads to: 



(r*)~ 3 / 2 (r*) a 
q 2( T *y 



1 



-1/2 



l-a 
v 3/2-2a 



(E10) 



This is result is inconsistent with the initial hypothesis 
M q ~ \j(p = l/q 2a+1 , unless a = 1/4, which is in con- 
tradiction with the initial hypothesis a < 1/4. Therefore, 
no values of a > 1/4 are authorized by the theory. 
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Then, we investigate the case 7 = 3/2 (or a = 1/4). 
Noting that Y 2 does not diverge for small times, we 
evaluate M q simply by taking the short time limit of the 
integrands in Eq. (E4): 



M a 



7r 2 q 2 J °° dr (^ v /r) _ 2e _ 



q 3/2> 

(Ell) 



where the last equali ty r esults from explicit integration 
and from the use of (E2). This result is consistent with 
the initial hypothesis M q = M^/q 3 ^ 2 . The exponent 
7 = 3/2 is the smallest that is authorized in the theory, 
and we therefore expect that the spectrum of the average 
reactive conformations decays as M q ~ 1/q 3 / 2 . 



Appendix F: Equations of the non-Markovian theory 

in 3D 



In this section, we show how to derive the equation ( 86 ) 



of the non-Markovian theory in the case of a space with 3 
dimensions. The quantities that are to be determined are 
the means m\ of the modes ai in the direction defined 
by the entrance angle of the reactive monomer to the 
target. Note that here we assume that the covariance 
of the modes at the reaction is equal to the stationary 
covariance, which in particular assumes the isotropy of 
the covariance matrix. More complicated equations could 
be derived if this assumption is released. 

Our starting relation is the following general integral 
equation, which is a reinterpretation of equation ( 83 ) and 
is the 3d equivalent of Eq. ( 26 ) : 



TP stat (R / )P stat (|a)|R / ) = 

J dt J ^[P(R / ,t|^,0)P(|a),t|R / ,t;^,0) 

- P(R /? t|{stat, P u r }, 0)P(|a), t|R/, t; {stat, P u r }, 0)] 

(Fl) 

We can assume that R/ = RfU z with < Rf < a with- 
out loss of generality. The self-consistent equations that 
define the moments mj will be derived by multiplying 
this last equation by the vertical component ai z of the 
jth mo de and by integrating over the conformations |a). 
More precisely, we first evaluate the following integral: 



/ 



d\a)a ir P(\a},t\RfU z , t; ttq, 0) 



m, ; e 



(Rir-Rf cos 6) (F2) 



This evaluation follows from the application of the prop- 
agation and projection formulas in the radial direction, 
by noting that the projection of RfU z over u r is Rf cos 0. 
Similarly, noting that the projection of RfU z over uq is 



—Rf sin#, we evaluate the following integral: 

bi(l -e~ Xit )R f s'mt 



I 



d|a)a^P(|a),t|Z f u 2 ,t;7TQ,0) = 



(F3) 



Noting that aj Z = a^ r cos# — a^sin^ and using the two 
evaluations ( F2|F3 ), we get: 



d\SL)Oi z P(\SL),t\RfU z ,t', TTq , 0) 



cos Omje 



TT-Xit 



(cosOR^-Rf) (F4) 
= cos 0^' u + R f bi(l - e- A * £ )/(A^) (F5) 



The same reasoning leads to the following relation for 
i > 2 (for which m ^ stat ' Xo > = 0): 



/ 



d\aL)a iz P(\sL),t\Rf\i z ,t] {Pou r , stat}, 0) = 



6i(l 



-Xit 



(cosORo-Rf) (F6) 



(FT) 



The explicit expression of P(RfU z , t\nn, 0) is: 

P(R f u z ,t\irn,0) = 

1 f (RfUz ~ R^r) 2 1 

[2^3/2 ex P| 2iP J 

We develop this expression at first order in Rf-. 

P(RfU z ,t\na,0) ~ 

-Rl/W) / R Rf 

' l + {hIlL cos + 0(R})\ (F8) 



[2tt^] 3 / 2 
Similarly, we have: 

P(R f \i z ,t\{&t&t, i? o u r },0) : 

e -(fio) 2 /(2>» RoRf 



4> 



[2tt^] 3 / 2 



cos8 + 0(R 2 f ) 



(F9) 



At this stage, we can multiply Eq. (JFTj) by a^, i ntegrat e 
over | a)) by using the intermediate expressions QF5|F6p . 
The next step consists in developpin g the r esult at first 
order in Rf by using the expressions flF8|F9| ). The lowest 
order term is proportional to cos and vanishes after the 
average over the angle 0. The term proportional to Rf 
is: 



r°° dt r 

Jo W 2 h 



dO sin x 



/i^P^cosfl) 



2 , hjl-e-^) 



e 2 ^ 



1 - (cos/9) — - — — -e 2 ^ 

w J ^ 



(F10) 



The last step of the calculation is a simple integration 
over 0. The result is exactly Eq. (86). 
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